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SUMMARY 


The  present  project  deals  with  the  fuel  spray  modelling 
in  multifuel  engines  using  the  droplet  behaviour  approach. 

The  multifuel  engine  is  able  to  burn  different  types  of  fuel 
having  various  physical  and  chemical  properties,  either  clear 
,  '  or  mixture  with  arbitrary  mixing  ratio.  Analysis  of  the  phy¬ 

sical  and  chemical  properties  affecting  mixture  creation  and 
f  combustion  have  been  carried  out  to  define  the  parameters  af¬ 

fecting  the  accuracy  of  spray  modelling.  It  is  found  that  the 
r  spray  is  affected  by  air  flow  field,  heat  transfer,  fuel  ato¬ 

mization,  droplet  behaviour,  fuel  evaporation  and  chemical 
kinetics.  Solution  of  the  flow  field  and  experimental  re- 
'  suits  of  heat  transfer  have  been  obtained. 

Experiments  were  performed  in  a  set  up  that  allowed 
variation  of  the  swirl  chamber  geometry  and  operating  condi- 
|  tions  to  evaluate  the  instantaneous  surface  heat  transfer  un¬ 

der  firing  conditions.  The  instantaneous  heat  flux  was  deter¬ 
mined  taking  into  consideration  the  instantaneous  composition 
of  products  and  variation  of  specific  heats.  A  generalized 
expression  has  been  established  for  the  calculation  of  heat 
transfer  to  the  chamber  wall  under  variable  design  and  operat¬ 
ing  conditions.^ 

Prediction  of  the  air  flow  field  is  based  on  the  solution 
of  the  finite  difference  analogue  of  the  governing  differen¬ 
tial  equations  for  the  transport  of  mass,  momentum  and  energy. 
The  k-e  model  of  turbulence  has  been  used.  A  special  tech¬ 
nique  has  been  established  to  satisfy  compressibility,  local 
continuity  and  reduce  the  mass  source  at  each  grid  node  to 
zero.  Measurements  of  air  flow  field  are  carried  using  a  hot 
wire  anemometer.  Agreement  between  predictions  and  measure¬ 
ments  is  fairly  good,  which  confirms  the  accuracy  of  the  pre¬ 
diction  procedure. 
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CHAPTER  L 


INTRODUCTION 


The  multifuel  engine,  which  is  able  to  burn  different 
types  of  fuel  either  clear  or  mixtures  with  arbitrary  mixing 
ratio,  is  required  to  ease  the  problem  of  disproportion  bet¬ 
ween  consumption  and  production  of  diesel  fuel.  It  must  be 
able  to  satisfy  burning  of  liquid  fuels  of  wide  distillation 
range,  and  preserve  the  same  power  output  and  specific  fuel 
consumption  at  various  speeds  and  loads.  Therefore,  the  fuel 
system  must  secure  reliable  delivary  of  fuel  into  engine,  good 
mixture  creation,  good  ignition,  smooth  burning  with  minimum 
emissions  and  minimum  corrosion  of  engine  elements.  In  addi¬ 
tion  easy  starting  of  engine  must  be  secured. 

The  multifuel  engine  has  been  known  at  the  beginning  of 
the  development  of  combustion  and  represented  by  the  engine 
IANZ-BULDOG  having  a  hot  bulb.  However,  some  design  difficul¬ 
ties  at  that  time  have  stopped  their  development  for  many 
years.  A  lack  of  some  kinds  of  fuel  for  combustion  during 
World  War  II  and  the  increase  of  diesel  fuel  consumption  due 
to  its  use  in  various  fields  have  urged  the  resumption  of  the 
development  of  multifuel  engine.  A  multifuel  piston  combustion 
may  be  designed  on  the  basis  of  either  a  gasolene  engine  or 
diesel  engine.  Howevc’-,  one  has  to  say  that  the  diesel  engine 
suits  more  the  multifuel  operation  since  it  works  with  higher 
compression  ratio  as  well  as  with  fuel  injection  in  the  cylin¬ 
der  which  appears  to  be  necessary  conditions  for  combustion  of 
heavier  kinds  of  fuel.  Moreover,  the  reliability  and  economy 
of  diesel  engines  have  been  acknowledged. 

In  the  early  stages  of  the  development  of  modern  high 
speed  diesel  engines  the  varieties  of  combustion  chamber  lay¬ 
outs  were  numerous.  There  is  no  doubt  that  this  early  inten¬ 
sive  work  resulted  in  better  stabilization  of  combustion  using 
two  types  of  combustion  chambers,  namely  the  open  chambers  a 
and  the  divided  ones  with  their  basic  three  categories;  the 
prechambers,  the  air  cells  and  the  swirl  combustion  chambers. 

The  multifuel  engine  must  be  prepared  with  a  combustion  chamber 
controlling  combustion  of  fuel  during  the  period  of  rapid  pres¬ 
sure  rise  to  obtain  better  burning  characteristics.  That  is 
why  most  of  designers  used  either  precombustion  chamber  or 
swirl  one  with  high  temperature  steel  inserts  to  preserve  the 
chamber  temperature  and  consequently  shorten  the  ignition  de¬ 
lay,  Actually  swirl  chambers  with  high  temperature  steel  in¬ 
serts  are  suitable  for  burning  multifuel. 

This  shows  why  the  air-fuel  mixing  process  prior  to 
autoignition  is  important  in  multifuel  engines  in  general  and 
in  swirl  chamber  multifuel  engines  in  particular.  The  principle 
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of  burning  in  swirl  chamber  being  injection  of  fuel  across  a 
rapidly  rotating  mass  of  air  in  the  combustion  chamber  in  such 
a  manner  that  the  air  has  to  find  the  fuel.  In  this  case  the 
mixing  process  includes  transport  of  fuel  from  the  spray  by 
the  swirling  air  flow,  deflection  of  the  fuel  spray,  impinge¬ 
ment  of  the  spray  on  the' combustion  chamber  wall  and  the  for¬ 
mation  of  turbulence  at  the  separation  boundaries  of  the  air 
Jet  in  the  swirl  chamber.  Besides,  the impinged  Jet  spread  on 
the  combustion  chamber  wall  forming  a  fuel  film  of  very  small 
thickness.  This  film  takes  its  required  heat  of  evaporation 
from  the  wall  and  surroundings  by  convection. 

The  analytical  solution  of  the  mixing  process  resulting 
from  the  space  and  film  mixture  creation  incide  swirl  combus¬ 
tion  chambers  is  complex.  Many  investigators  have  studied 
the  problem  of  mixture  formation  and  combustion  in  diesel  en¬ 
gines,  but  most  of  them  treated  the  problem  either  in  general 
or  by  formulating  some  emperical  relations  to  represent  the 
individual  processes.  Only  few  investigators  have  tried  to 
formulate  the  mixture  formation  process  and  combustion  in  the 
swirl  chamber  on  the  basis  of  some  known  physical  laws. 

Along  this  line  the  present  project  is  conducted  to  deal 
with  a  theoretical  model  and  experimental  verification  of  the 
fuel  spray  trajectory  in  multifuel  engines  using  swirl  combus¬ 
tion  chambers.  This  spray  modelling  is  required  to  investigate 
speedily  and  with  reasonable  details  and  accuracy  the  combus¬ 
tion  of  fuel  having  various  physical  and  chemical  properties 
and  to  study  the  development  of  multifuel  engines.  The  fuel 
mixture  creation  and  combustion  are  generally  affected  by  the 
air  flow  field, heat  transf er, fuel  atomization,  droplet  beha¬ 
viour,  fuel  evaporation, chemical  kinetics,  delay  period  and 
combustion  chamber  shape.  The  mathematical  model,  which  in¬ 
cludes  all  these  factors,  can  give  quite  sufficient  accurate 
results.  These  factors  have  been  solved  in  successive  stages. 

In  this  first  Annual  Technical  report,  analysis  of  the 
physical  properties  affecting  mixture  creation  and  combustion 
in  multifuel  engines,  solution  of  the  air  flow  field  and  heat 
transfer  in  swirl  chambers  have  been  dealt  with.  The  govern¬ 
ing  conservation  equations  of  the  dimensional  flow  field  in 
the  swirl  chamber  have  been  formulated.  The  k-e  model  of 
turbulence  has  been  used  to  define  the  turbulent  viscosity  and 
the  turbulent  exchange  coefficient.  The  equations  are  solved 
using  the  finite  difference  technique  in  which  the  compres¬ 
sibility  and  local  continuity  were  satisfied.  A  special  tech¬ 
nique  has  been  established  to  reduce  the  mass  source  at  each 
grid  node  to  zero.  An  experimental  plant  has  been  built  to 
measure  the  flow  field  in  an  actual  swirl  chamber  using  the 
hot  wire  anemometer.  A  comparison  over  a  wide  range  of  data 
is  required  to  validate  the  mathematical  model.  The  solution 


of  the  flow  field  has  been  conducted  by  Prof.  ELKOTB, M.M.  , 

Dr.  ABOU-ALLEIL ,  M.M.  and  Eng.  SALEM,  S. 

The  instantaneous  surface  heat  transfer  in  swirl  combustion 
chamber  of  various  geometry  has  been  determined  experimentally 
under  firing  conditions.  Generalized  expressions  have  been 
established  for  the  calculation  of  heat  transfer  to  the  cham¬ 
ber  wall  under  variable  design  and  operating  conditions. 

Prof.  ELKOTB ,  M.M. ,  Dr.  AWAD,  M.M.  and  Eng.  ESMAIL  are  res¬ 
ponsible  for  this  part  of  the  report. 

■  •  The  analysis  of  the  influence  of  fuel  properties  on  mix¬ 
ture  creation  and  consequently  the  engine  performance  together 
with,  the  prelimenary  measurement  of  engine  performance  when 
burning  multifuel  have  been  conducted  by  Prof .ELKOTB,  M.M.  and 
Dr.  ELGOHEIMY ,  M.N. 

Experimental  set  ups  for  the  investigation  of  other  men¬ 
tioned  important  parameters  have  been  built  up.  Their  details 
together  with  the  results  will  be  discussed  in  other  Technical 
Reports . 


CHAPTER  2 


INFLUENCE  OF  PHYSICAL  PROPERTIES  OF  FUEL  ON 
MIXING  AND  COMBUSTION  PROCESSES 


2-1  INTRODUCTION 


The  multifuel  engine  must  satisfy  the  common  demands  of 
all  engines  such  as  high  specific  power  output,  high  economy, 
minimum  size  and  weight,  high  reliability  and  easy  maintenance. 
Moreover,  it  must  be  able  to  meet  further  important  demands  as 
follows  : 

a-  to  enable  combustion  of  liquid  fuel  of  both  natural  and 
synthetic  origin  of  wide  distillation  range,  +  30°C  to  450°C. 
i.e  gasoline.  Kerosene,  aviation  fuel,  Diesel  fuel  and  even 
all  kinds  of  lubricating  oils. 

b-  to  preserve  approximately  the  same  power  output  and 
operating  characteristics  for  all  fuels. 

At  the  same  time  the  multifuel  engine  must  operate  in  the  same 
way  when  using  either  clear  fuels  or  mixtures  with  any  arbit¬ 
rary  ratio.  Table  2.1  shows  the  physical  and  chemical  proper¬ 
ties  of  main  Fuels  produced  in  Egypt.  It  is  evident  from  this 
table  that  the  fuels  to  be  used  in  multifuel  engines  posses 
various  physical  properties  which  affect  the  mixing  and  com¬ 
bustion  processes  differently.  Therefore,  the  fuel  system  must 
secure  reliable  aelivary  of  fuel  into  engine,  good  mixture 
creation,  good  ignition,  smooth  burning  with  minimum  emission 
and  minimum  corrosion  of  engine  elements.  In  addition  easy 
starting  of  engine  must  be  secured.  These  demands  being  laid 
upon  multifuel  are  not  small  and  their  successful  fulfuilment 
strikes  at  present  with  several  difficulties.  A  brief  analysis 
of  various  physical  properties  of  the  mentioned  fuel  will  be 
carried  out  to  explore  their  effects  on  the  engine  performance. 


2-2  ANALYSIS  OF  VARIOUS  PHYSICAL  PROPERTIES 

a-  Low  distillation  temperature  of  light  fuels  slight  fuel 
such  as  gasolene  and  aviation  fuel  will  cause  vapour  lock  due 
to  evaporation  of  light  fractions.  This  may  happen  in  the 
low  pressure  branch  of  the  fuel  system  when  exposed  to  high 
temperatures  from  the  engine  and  causes  an  irregular,  inter¬ 
rupted  and  incorrect  metering  of  fuel.  Restarting  of  warm  en¬ 
gine  will  be  interrupted  when  using  light  fuels  as  the  fuel 
will  evaporate  as  a  result  of  the  high  temperatures  achieved. 
The  temperatures  of 'the  fuel  system  may  reach  a  value  of  about 
50°C  or  more  (  up  to  60°C)  exceeding  the  initial  distillation 
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temperat-  .  o.  This  difficulty  of  operation  increases  with  the 
increase  of  the  atmospheric  temperature,  the  condition  of  ope¬ 
ration  most  prevailing  in  Egypt. 

The  temperature  of  the  fuel  delivery  parts  must  not  ex¬ 
ceed  40  °C  if  vapour  lock  is  to  be  prevented.  More  effective 
fuel  delivery  pump  which  delivers  greater  quantity  of  fuel, 
is  required  in  multifuel  engines  as  the  flow  of  fuel  cools  the 
fuel  pump  intensively.  Thus, to  secure  the  gooc’  ability  of 
restarting  of  warm  engine  it  would  be  benifitial  to  equip  the 
fuel  system  with  an  electric  delivery  pump,  which  is  independent 
on  the  drive  from  the  engine.  The  increase  of  the  delivery 
pressure  of  fuel  to  1.5  up  to  3  bar  will  also  limit  the  tenden¬ 
cy  to  create  vapour  lock. 

b-  Low  viscosity  of  lighter  fuels. 

Low  viscosity  of  fuel  is  obtained  when  using  light  fuels. 
This  causes  leakage  of  fuel  through  clearances  in  injection 
units  resulting  in  incorrect  metering  of  fuel.  In  addition  ex¬ 
cessive  atomization  of  fuel  is  obtained  causing  small  fuel  penet¬ 
ration,  lack  of  homogenity  and  burning  in  the  vicinity  of  the 
injector.  This  may  cause  overheating,  carbonization  and  de¬ 
formation  of  some  parts.  Moreover,  the  lubrication  of  the  mov¬ 
ing  elements  of  the  pump  decreases  causing  excessive  wear. 

Leakage  of  light  fuels  can  be  prevented  either  by  oil 
sealing  rings  or  by  manufacturing  special  grooves  on  plunger  in 
which  the  fuel  may  sedimentate  and  then  gets  drained  into  the 
suction  part  of  the  pump.  The  injection  pump  must  also  be  supp¬ 
lied  .with  lubricating  oil  from  the  engine  to  decrease  the  wear 
of  i  s  moving  parts. 

c-  Low  value  of  cetane  number. 

Light  fuels  have  low  values  of  cetane  number  in  compari¬ 
son  with  diesel  fuel,  creating  the  following  difficulties  : 

i,  high  self  ignition  of  fuel  as  well  as  high  delay  period. 

ii,  hard  running  and  difficult  starting  of  the  engine  due 
to  small  compression  temperature. 

The  cetane  number  of  the  fuel  depends  on  both  the  chemi¬ 
cal  and  fractional  composition  of  fuel.  The  alkaline  hydro¬ 
carbons  have  the  highest  cetane  number  while  the  aromatie  hyd¬ 
rocarbons  have  the  lowest  cetane  number.  Cyclones  have  middle 
values  between  alkalines  and  aromatics.  Lighter  fuels  have 
small  cetane  number  resulting  in  difficult  ignition  of  fuel 
when  used  in  diesel  engines.  However,  high  cetane  fuels  posses 
low  ignition  delay  and  burn  in  the  vicinity  of  the  injector. 

That  is  why  the  air  present  far  away  from  the  injector  does  not 
take  part  in  the  burning  process  and  consequently  incomplete 
burning  of  fuels  may  occur. 
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This  causes  smoky  combustion  and  drop  of  engine  economy. 

The  problems  of  hard  running  and  long  delay  period  are 
one  of  the  main  difficulties  of  multifuel  engines  working  with 
various  liquid  fuels.  The  analysis  of  these  problems  and  the 
possibility  of  their  solution  will  be  dealt  with  later  on. 

c-  Saturation  fuel  pressure. 

It  is  known  that  the  saturation  pressure  of  the  fuel 
characterizes  the  evaporation  of  fuel.  The  REID  pressure  of 
the  diesel  fuel  is  zero  while  it  may  achieve  a  value  of  100 
torrs  for  petrol  fuel.  Thus  the  saturation  pressure  increases 
with  the  increase  of  light  fractions.  The  increase  of  satura¬ 
tion  pressure  causes  an  excessive  creation  of  vapour  lock. 
Moreover,  small  amount  of  heat  is  required  to  evaporate  the 
same  quantity  of  fuel,  it  is  required  to  use  fuels  of  REID 
pressure  not  exceeding  450  torr. 

d-  fractional  composition  of  fuels. 

The  fractional  composition  of  fuels  affects  directly  the 
ability  of  fuel  to  evaporate, atomize  and  disperse.  The  time 
of  evaporation  in  high  speed  diesel  engines  is  very  short. 

It  is,  therefore,  necessary  for  both  light  and  heavy  fractions 
contained  in  diesel  engine  fuels  to  be  in  a  suitable  ratio. 
Heavy  components  of  fuel  evaporate  slowly  and  imperfect  com¬ 
bustion  is  obtained  leading  to  smoke  and  soot  concentration  in 
the  exhaust  gases.  Consequently,  the  engine  efficiency  dec¬ 
reases. 

While  .evaporation  of  light  components  is  quicker,  their  point 
of  inflammation  is  higher  than  that  of  heavier  components.  A 
greater  quantity  of  light  components  in  fuels  may  cause  high 
fuel  concentration  in  the  zones  behined  the  injector.  This 
may  cause  high  increase  of  pressure  rate  and  overheating,  re¬ 
sulting  in  hard  running  of  engine. 

e-  change  of  calorific  value  and  specific  weight  of  various 

fuels . 

The  change  of  calorific  value  is  not  substantion  as  it 
does  not  exceed  3%.  However,  the  specific  weight  varies  by 
10  up  to  15%.  The  decrease  of  the  specific  weight  causes  a 
decrease  of  ihe  quantity  of  fuel  injected  and  consequently  a 
decrease  of  the  output  power  from  the  engine  and  engine  economy. 

Therefore,  the  injection  pump  of  multifuel  engine  must 
be  equipped  with  an  arrangement  for  changing  the  injected  quan¬ 
tity  independence  of  the  fuel  specific  weight. 

The  effect  of  fuel  density  change  is  considerably  lower 
than  that  of  viscosity.  However,  they  must  be  considered  to¬ 
gether  with  their  change  in  the  operating  temperature  when 
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selecting  the  holes  of  the  injectors. 

2-3  MULTIFUEL  COMBUSTION  IN  DIESEL  ENGINES 


Combustion  of  fuels  in  diesel  engines  is  a  very  compli¬ 
cated  problem.  Fuel  inflammation  occurs  from  the  heated  up 
air  in  the  combustion  chamber.  For  that  reason  oompression 
ratio  of  more  than  13  is  necessary  to  create  convenient  con¬ 
ditions  for  burning  before  TDC.  The  fuel  is  injected  in  the 
combustion  chamber  with  an  injection  advance  ranging  between 
20  to  40  degrees  before  TDC  according  to  the  engine  speed. 
Thereafter  the  ignition  takes  place  after  a  period  of  time 
equivalent  to  10-20  degrees  of  the  crankshaft,.  The  variation 
of  gas  pressure  inside  the  cylinder  with  crankangle  for  gaso¬ 
lene  and  diesel  fuel  is  shown  in  Fig. (1) .  The  preparation  of 
fuel  for  combustion  is  responsible  for  the  ignition  delay  pe¬ 
riod,  The  whole  process  of  burning  as  well  as  the  performance 
of  the  engine  depends  on  the  ignition  delay  period.  The  quan-r 
tity  of  fuel  delivered  during  the  delay  period  wiibl  start  to 
burn  directly  after  creation  of  active  centres.  It  follows 
that  the  longer  the  delay  period  is  the  more  quantity  of  fuel 
is  accumulated  in  the  combustion  space  and  then  rapid  increase 
of  pressure  is  obtained.  Hard  operation  of  the  engine  is  ob¬ 
tained  when  the  rate  of  pressure  rise  (roughness  of  combustion) 
exceeds  6  bar/degree  crank. 

This  increase  of  pressure  rise  rate  is  accompanied  by  exces¬ 
sive  wear  of  bearing,  crankshaft  and  rumbling  sound.  High 
harpl  operation  of  the  engine  is  also  occompanied  by  smoke  out¬ 
let,  and  high  specific  fuel  consumption. 

From  the  above  can  be  seen  that  increasing  the  delay 
period  increases  the  hardness  of  diesel  engine  work  while  dec¬ 
reasing  it  makes  the  engine  operation  softer.  However,  hav¬ 
ing  too  small,  a  delay  period  and  at  the  same  time  long  fuel 
delivery  period,  lead  to  imperfect  combustion,  and  soot  is 
obtained  in  the  exhaust.  The  ignition  delay  period  is  deviu- 
ed  into  physical  and  chemical  delay  periods  which  overlap  each 
other , 


The  fuel  atomization,  partial  evaporation  and  heating 
up  of  vapours  to  the  inflammation  temperature  occur  during 
the  physical  delay  period.  The  formation  of  active  centres 
and  oxidation  occur  during  the  chemical  delay  period.  The 
physical  ignition  lag  has  a  considerable  value  when  using 
heavy  fuels.  It  is  evident  from  Fig.l  that  the  ignition  de¬ 
lay  of  gasolerye  burnt  in  diesel  engine  is  considerably  higher 
than  that  of  diesel  fuel  and  consequently  hard  operation  is 
obtained.  Therefore,  a  normal  diesel  engine  is  not  suitable 
to  work  With  gasolene.  The  multifuel  engine  must  satisfy  the 


injection  duration 


FIG  .1 


BEHAVIOUR  OF  DIESEL  ENGINE  WITH  DIESEL  FUEL 
AND  GASOLENE. 
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same  pressure  variation  as  shown  in  Pig. 2.  A  method  must  be 
find  to  remove  the  quick  increase  of  pressure  of  the  uncontroll¬ 
ed  part  (c'd')  on  the  diagram.  Therefore,  it  is  required  to 
control  the  ignition  delay  period.  To  try  to  solve  this  prob¬ 
lem  it  is  required  to  define  accurately  the  influence  of  va¬ 
rious  factors  on  the  delay  period. 

The  delay  period  is  influenced  by  the  kind  of  fuel,  oxy¬ 
gen  concentration,  concentration  of  residual  gases  and  presence 
of  catalyzers  and  admixtures. 

The  parafinic  hydrocarbons  have  the  shortest  ignition  de¬ 
lay  due  to  the  long  chain  structure  while  aromatics  have  the 
longest  ignitic  1  delay  due  to  their  ring  structure.  Therefore, 
in  comparison  with  the  standard  diesel  fuel  the  ignition  of 
Kerosene  is  somewhat  greater,  its  increase  of  pressure  smooth¬ 
er  and  maximum  pressure  smaller.  On  the  contrary,  gasolene 
(aromatic)  has  a  long  delay  period  thus  increase  of  medium  pres¬ 
sure  is  sharp  and  the  maximum  pressure  is  high.  The  self-igni-. 
tion  temperature  of  various  kinds  of  fuel  were  found  bo  be 

gasolene  470  ~  530  nC 

Jet  fuel  (aviation  kerosene)  240  -  380  °C 

diesel  fuel  250  -  270  °C 

Fuels  with  low  self  ignition  temperatures  such  as  kerosene 
and  diesel  fuel  are  suitable  for  burning  in  diesel  engines. 

The  increase  of  oxygen  concentration  in  mixture  from  20% 
to  33%  decreases  the  ignition  delay  from  80xl0~  sec  to  55x10  4 
sec.  The  change  of  ignition  delay  with  oxygen  concentration  is 
approximately  linear.  Therefore  the  use  of  oxygen  improves  the 
engine  starting  as  well  as  combustion. 

It  is  also  known  that  a  great  number  of  physical  factors  * 
are  affecting  the  delay  period  such  as  evaporation  characteris¬ 
tics  of  fuel,  its  viscosity,  air  temperature  at  the  end  of 
compression,  fresh  charge  temperature  and  density,  air  swirl 
and  heat  transfer  in  combustion  chamber.  For  high  speed  en¬ 
gines  heavy  fuel  has  no  time  to  evaporata. 

The  increase  of  turbulence  intensity  increases  the  quan¬ 
tity  of  heat  transferred  to  the  cylinder  wall,  resulting  in  a 
decrease  of  pressure  and  temperature  of  gases  at  the  end  of 
compression  stroke.  This  decrease  of  temperature  increases 
the  delay  period.  At  the  same  time  the  increase  of  turbulence 
intensity  improves  the  fuel  atomization,  heat  transfer  to  fuel 
droplets  and  evaporation,  resulting  in  a  decrease  of  the  delay 
period.  *  Therefore,  it  is  possible  to  say  that  the  ignition 
delay  is  shortened  with  the  increase  of  the  turbulence  inten¬ 
sity. 


110  bar 


COMPARISON  BETWEEN  DIESEL  ENGINE  AND  MULTIFUEL 
ENGINE  WORKING  WITH  DIESEL  FUEL  AND  GASOLENE, 
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The  injection  pressure  is  one  of  the  factors  that  deter¬ 
mines  the  quality  of  fuel  atomization.  The  increase  of  fuel 
pressure  decreases  the  fuel  droplets  diameter  and  consequent¬ 
ly  the  total  surface  of  fuel  drops  increases.  This  increases 
the  heat  exchange  to  fuel  and  relative  oxygen  concentration 
on  the  surface  of  fuel  droplets,  thus  accelerating  the 
physical-chemical  preparation  of  fuel  for  ignition.  The  law 
of  fuel  delivery  has  no  effect  on  the  ignition  delay.  The 
change  of  fuel  delivery  causes  a  change  of  the  fuel  delivered 
per  crank  angle  and  leads  to  the  change  of  heat  release  rate 
and,  consequently,  the  pressure-time  diagram.  The  number  and 
size  of  injecting  nozzle  as  well  as  the  location  of  the  in¬ 
jector  have  a  substantial  effect  on  the  fineness  and  unifor¬ 
mity  of  atomization.  However ,  there  will  be  a  substantial 
effect  on  the  spray  length  inside  the  combustion  space. 

Finally  it  is  evident  that  gasolene  burning  in  diesel 
engines  requires  an  increase  of  compression  ratio,  higher  ope¬ 
rating  temperatures,  avoidance  of  creation  of  vapour  lock  and 
regulation  of  fuel  quantity. 


2-4  CONTROLLING  COMBUSTION  OF  MULTIFUEL 

From  the  analysis  of  the  factors  affecting  the  delay  pe¬ 
riod  it  is  evident  that  the  compression  ratio,  combustion 
chamber  shape  and  design,  method  of  mixture  creation  and  init¬ 
ial  conditions  to  engine  have  valuable  importance  when  design¬ 
ing  the  multifuel  engine.  Moreover,  the  control  of  the  course 
of  'burning  has  a  direct  effect  on  the  performance  of  multifuel 
engine.  Controlling  the  course  of  burning  to  secure  multifuel 
combustion  has  been  satisfied  by  shortening  the  delay  period, 
regulating  the  combustion  of  fuel  during  the  rapid  pressure 
rise  period  and  mixture  creation  in  the  combustion  chamber  space. 
That  is  why  most  of  the  designers  used  either  precombustion 
chambers  as  shown  in  Fig. 3  or  swirl  chambers  as  shown  in  Fig. 4. 
These  chambers  are  equipped  with  hot  surfaces  and  thermal  in¬ 
sulations  to  preserve  the  chamber  temperature  at  high  values 
sufficient  to  shorte.i  ignition  delay  of  armmatic  fuels.  Pre¬ 
combustion  chambers  are  used  by  Daimler-Benz  and  M.A.N  fac¬ 
tories.  Swirl  combustion  chambers  are  used  by  Deutz  in  their 
multifuel  engines.  The  English  multifuel  eggines  are  in  most 
cases  two  stroke  engines  with  opposite  pistons  (Leyland,  Rolls- 
Royce,  commer  and  country  climax) . 

Air  cell  wi  two  holes  injector  can  also  be  used.  Paet 
of  the  fuel  is  directed  tangentially  on  the  wall  while  the 
i  other  .part  is  directed  towards  the  centre.  Such  combustion 

chamber  is  used  by  MAN  Factory,  Fig. 5.  Five  percent  of  the 
injected  fuel  is  directed  tangentially  on  the  wall  of  the  sphe¬ 
rical  combustion  chamber  located  in  the  piston. 


injector 


spark  plug 


a)  external  channel  of  casing 

b)  internal  channel  of  casing 


fig. 3 


PRECOMBUSTION  CHAMBER  (MWM)  FOR  MULTIFUEL  ENGINE 
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2-5  PRELIMENARY  RESULTS  OF  BURNING  VARIOUS  FUELS  IN  A 

DIESEL  ENGINE 

It  was  clear  from  the  previous  analysis  that  the  engine 
performance  was  influenced  by  various  factors  affecting  the 
mixture  creation  during  the  delay  period.  It  is  therefore  re¬ 
quired  to  measure  the  engine  performance  when  burning  either 
clear  fuel  or  mixtures  of  arbitrary  ratio.  This  is  required 
to  define  the  object  of  this  work. 

An  Egyptian  made  diesel  engine,  Nasr  92-0501,  was  used 
for  this  purpose.  It  is  a  stationary,  water  cooled,  four  stroke, 
2-cylinder  engine.  The  main  data  of  this  engine  are  given  in 
Table  2.  The  cylinder  head  is  provided  with  a  precombustion 
chamber.  Pintle  type  injector  with  opening  pressure  of  130 
bar  was  used  to  inject  the  fuel  at  an  injection  advance  angle 
of  18  degrees  before  TDC.  The  engine  was  equipped  with  the 
required  instrument  for  measuring  the  operational  condition 
and  the  engine  performance.  A  hydraulic  dynamometer  was  used 
to  measure  the  brake  power.  The  exhaust  gas  temperature  was 
measured  by  a  Fe-constantan  thermocouple.  The  measurement  of 
smoke  of  exhaust  gases  was  carried  out  by  the  through-flow  gas 
analyser  type  Hartrige. 

The  tests  have  been  carried  out  with  the  following  fuels; 
diesel  fuel  p  =  0.842  gm/cm  ,  gasolene  70  O.N  p  =  0.732  gm/cm. 
Kerosene  p  =  0.8115  gm/cm3,  jet  Kerosene  p  =  0.811  gm/cm3  and 
mixture  of  diesel  oil  and  gasolene  p  =  0.789  gm/cm3.  The  in¬ 
jection  pump  was  adjusted  according  to  the  producer  recommenda¬ 
tion  to  give  8.5  cm3  for  100  strokes  at  a  speed  of  500  RPM. 

Some  of  the  results  obtained  for  the  performance  of  the  engine 
when  burning  the  mentioned  fuels  are  illustrated  in  Fig. 6.  It 
was  found  that  the  power  decreased,  the  specific  fuel  con¬ 
sumption  increased  and  the  smoke  decreased  with  the  increase  of 
the  gasolene  ratio  in  the  diesel  fuel.  It  was  also  found  that 
small  changes  occured  when  using  either  kerosene  or  jet  kerosene 
as  shown  in  Fig. 7. 

Table  2  Nasr  592-0501  Diesel  Engine  Data 

Type  Four  stroke,  precombustion  chamber,  water  cooled  natu¬ 
rally  aspirated  diesel  engine. 

Bore  and  stroke  125  mm  and  140  mm 
Capacity  3.44  litres 

Compression  ratio  22 

Max.  output  34  HP  at  1500  RPM 

Max.  torque  17  kpm 

Injection  advance  18  degrees  before  TDC. 


CHAPTER  3 


OBJECT  OF  THE  PRESENT  WORK 


It  is  clear  from  the  previous  review  that  the  modified 
diesel  engine,  which  uses  swirl  chamber,  is  suitable  to  bum 
smoothly  multifuel,  and  the  performance  of  this  engine  de¬ 
pends  mainly  on  mixture  preparation.  Many  researches  have 
been  carried  out  to  investigate  the  nature  of  combustion  in 
swirl  chamber  diesel  engine.  Some  of  these  studies  have  been 
carried  out  on  experimental  basis  to  determine  the  effect  of 
swirl  chamber  shape on  the  engine  perf ormance .  These  studies 
were  aimed  at  finding  out  the  swirl  chamber  dimensions  for 
optimum  mixture  formation  to  satisfy  minimum  fuel  consumption, 
low  roughness  of  combustion,  small  noise,  easy  starting  and 
high  outputs.'  Most  of  these  researches  led  to  a  rough  eva¬ 
luation  of  the  swirl  chamber  dimensions.  Only  few  of  these 
researches  {  1  -  6  {  have  tried  to  establish  a  convenient 

theoretical  model  of  the  space-film  mixture  creation.  How¬ 
ever,  this  model  is  affected  by  the  turbulence  motion  of  air, 
heat  transfer  to  the  chamber  and  droplets,  droplet  size  dis¬ 
tribution,  rate  of  fuel  evaporation  and  the  rate  of  chemical 
reaction.  Moreover,  the  experimental  validity  of  the  model 
has  not  been  yet  done. 

The  r  •  be  of  heat  release  and  the  roughness  of  combus- 
cion  are  controlled  by  the  rate  of  evaporation.  The -mixture 
formation  inside  the  swirl  chamber  depends  mainly  on  the  air 
velocity  pattern,  the  behaviour  of  the  fuel  spray  either  in 
space  or  on  wall  and  the  mechanism  of  chemical  reaction. 

The  air  velocity  pattern  was  found  to  depend  upon  the 
swirl  chamber  dimensions  and  the  pressures  in  the  main  and 
swirl  chambers.  The  pattern  of  air  movement  inside  the  swirl 
chamber  was  found  to  be  a  solid  vortex  motion  at  the  central 
part  and  free  vortex  motion  at  the  periphery  j  4  { .  This 
has  facilitated  the  establishment  of  a  satisfactory  theore¬ 
tical  model  for  the  mean  tangential  air  velocity  and  medium 
pressure  inside  the  sviirl  chamber  at  various  crank  angles. 
However,  the  determination  of  the  average  velocity  does  not 
reveal  the  details  of  interaction  between  fuel  and  air. 

Over  the  past  decade  a  considerable  effort  has  been  expended 
on  computational  procedure  with  the  aid  of  mathematical  models 
of  turbulence  for  the  determination  of  air  velocity  in  recip¬ 
rocating  engines  |7-9 { .  These  prediction  procedures  are  bas¬ 
ed  mainly  on  the  solution  of  the  finite  difference  of  the 
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governing  differential  equations  for  the  transport  of  mass 
momentum  and  energy.  However,  no  explicit  consideration  of 
the  physical  properties  in  swirl  chamber  has  been  taken  in 
these  attempts.  These  include  the  continuous  variation  of 
medium  pressure  and  temperature,  heat  transfer  during  flow 
and  the  shape  of  the  combustion  chamber.  Moreover,  no  ex¬ 
perimental  measurements  of  the  air  flow  velocity  were  car¬ 
ried  out  to  carefully  compare  their  values  with  the  predicted 
one  in  order  to  obtain  reliable  information  on  the  constants 
of  the  tubbulence  model. 

The  space  fuel-air  mixing  has  been  studied  by  several 
investigators  1 5,6,1].  ELKOTB  et  al  {  2  j  studied  the 
space  fuel-air  mixing.  Their  study  was  based  on  the  beha¬ 
viour  of  a  droplet  of  mean  diameter  moving  at  the  fuel  spray 
front  along  the  spray  axis  and  subjected  to  variable  heat 
transfer  and  drag  from  the  surrounding  medium.  Thereafter 
they  proposed  a  somewhat  miodif  ied  model  which  approximates 
the  fuel  spray  with  a  finite  number  of  droplet  size  groups. 

They  discussed  various  factors  affecting  the  spray  formation 
and  obtained  the  boundary  conditions  of  the  spray  at  the  in¬ 
stant  of  impingement.  Moreover,  ELKOTB  1101  established  a  theo¬ 
retical  model  for  the  fuel  film  motion  on  the  surface  of  a 
swirl  chamber.  But  his  work  is  merely  theoretical  and  need 
be  verified  experimentally  and  extended  to  cover  different 
factors  affecting  the  wall  jet. 

Data  on  drop  size  from  various  types  of  fuel  injectors 
and  data  on  single  droplet  vaporization  rate  are  plentiful. 

But  Only  limited  data  have  been  published  on  spatial  fuel  dis¬ 
tribution  and  vaporization  rates  of  sprays.  No  data  have 
bean  published  on  spatial  fuel  distribution  and  vaporization 
rates  of  mixed  fuel  with  various  ratios .  Such  information  is 
necessary  for  the  development  of  the  spray  model  in  swirl 
chamber  diesel  engines. 

From  the  previous  discussions,  one  can  conclude  that 
the  following  parameters  are  the  main  factors  controlling 
the  space-film  mixture  formation  in  swirl  chamber  diesel  en¬ 
gines. 

, 1-  The  operating  conditions  and  their  variations  with 
crank  angle,  swirl  chamber  geometry,  and  heat  transfer  to 
combustion  chamber  walls. 

2-  Air  velocity  pattern  inside  the  swirl  chamber  and  its 
change  with  different  factors. 

3-  The  mechanism  of  spray  formation  either  in  space  or  on 
the  wall  and  the  factors  affecting  it. 
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4-  The  mechanism  of  film  heating,  evaporation  and  mixing 
with  air  to  form  combustible  mixture  and  the  factors  con¬ 
trolling  these  processes. 

5-  The  spatial  droplet  size  distribution  of  multifuel  and 
factors  affecting  it. 

6-  Droplet  behaviour  and  its  evaporation  for  fuel  mixtures 
and -factors  affecting  it. 

7-  Velocity  of  chemical  reaction  and  chemical  kinetics  of 
fuel  mixtures  and  factors  affecting  it. 

The  present  project  intended  to  reveal  more  information 
about  the  modelling  of  multifuel  spray  formed  inside  swirl 
chamber  diesel  engines.  The  study  of  the  mentioned  parameters 
are  required  to  satisfy  a  sufficient  accuracy  of  the  spray 
modelling.  Several  experimental  setups  have  been  built. to 
investigate  these  factors. 

A  set  up,  in  which  the  swirl  chamber  geometry  can  be 
changed,  has  been  built  to  study  the  spatial  variation  of 
heat  transfer  in  swirl  chambers.  These  experimental  results 
will  be  used  for  the  validation  of  the  theoreticalrresults . 

Another  experimental  set  up  consisting  of  transparent 
swirl  combustion  chamber,  injection  system,  high  speed  photog¬ 
raphy,  hot  wire  anemometer  and  the  necessary  equipment  for  re¬ 
cording  various  parameters  has  been  built.  This  is  required 
for  the  measuring  of  air  velocity  and  spray  behaviour.  A 
third  one„has  been  built  to  study  the  spatial  atomization  of 
fuel  mixtures  using  the  slide  sampling  technique. 

Ohher  interesting  parameters  are  the  droplet  behaviour 
and  fuel  evaporation.  An  experimental  set  up  was  built  to 
study  the  Former  parameter  using  the  flight  droplet  technique. 
Another  setup  was  built  to  measure  the  latter  one  using  iso*- 
kinetic  sampling  probe  to  pick  up  the  evaporated  fuel,feurn  it 
and  analyse  it  by  infrared  analyzers.  This  can  give  the  rate 
of  fuel  evaporation. 

Moreover,  it  is  planned  to  study  the  chemical  kinetics 
of  the  multifuel,  and  a  shock  tube  has  been  built  for  this  pur¬ 
pose. 

In  the  present  report  it  is  intended  to  present  some 
new  results  for  the  heat  transfer  in  swirl  combustion  chambers. 
The  instantaneous  heat  flux  was  determined  and  a  generalized 
expression  has  been  established  for  the  calculation  of  the 
heat  transfer  to  the  swirl  chamber  wall.  Moreover,  the  pre¬ 
diction  and  measurements  of  flow  velocity  in  swirl  chambers 


are  presented.  Measurements  of  air  velocity  are  mainly  carried 
out  by  a  hot  wire  anemometer.  Prediction  of  air  velocity  are 
based  on  the  solution  of  the  finite  difference  form  of  the 
governing  differential  equations  for  the  transport  of  mass, 
momentum  and  energy. 
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CHAPTER  A 


SPATIAL  VARIATION  OF  HEAT  TRANSFER 
IN  SWIRL  COMBUSTION  CHAMBERS 


ABSTRACT 


Experiments  were  performed  to  evaluate  the  instantane¬ 
ous  surface  heat  transfer  in  a  diesel  engine  *wirl  chamber 
under  firing  conditions.  An  experimental  set  up,  in  which 
the  swirl  chamber  facilitated  the  variation  of  the  relative 
swirl  chamber  volume,  relative  connecting  channel  area,  in¬ 
clination  angle  of  the  connecting  channel  and  fuel  pressure 
had  been  built.  The  instantaneous  gas  temperature  was  deter¬ 
mined  using  the  determined  coarse  of  heat  release,  coarse  of 
pressure  in  the  swirl  chamber  and  the  equation  of  state. 
Thereafter,  the  instantaneous  heat  flux  was  determined  taking- 
into  consideration  the  momentary  composition  of  products  and 
variation  of  specific  heat.  Measuring  the  mean  wall  tempera¬ 
ture  together  with  the  previous  mentioned  results  it  was 
possible  to  determine  the  instantaneous  heat  transfer  coeffi¬ 
cient. 

Generalized  expression  has  been  established  for  the 
calculation  of  the  heat  transfer  to  the  swirl  chamber  wall 
under  variable  design  and  operating  conditions. 


NOMENCLATURE 

a-^-a^  constants  of  eqn.  (3) 

ag  relative  swirl  chamber  volume,  Vs/Vc,  dimensionless 
Ae  area  of  tangential  port,  m^ 

As  area  of  swirl  chamber  surface,  nt2 
Ca  discharge  coefficient,  dimensionless 

Cmp  mean  piston  speed,  m/s 

Cp  specific  heat  of  gas  at  const,  pressure,  J/Kg  K 
d  diameter,  m 

D  engine  bore  diameter,  m 

Eu  Euler  Number,  Pf0/p  C2 

g  gravity  acceleration,  m/s2 

h  convective  heat  transfer  coefficient,  J/m2h  K 

K  thermal  conductivity  coefficient,  J/m  h  K 

m  mass  of  fuel,  Kg/cycle. 

M  mass  of  gas,  Kg 

n  poly tropic  index,  dimensionless 

N  engine  speed,  RPM 

Nu  Nuselt  Number,  hD/K 

P  pressure,  bat. 
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Pr  Prandtl  Number,  Cp  y/K 

Qr  heat  release,  KJ 

?w  heat  transfer  to  wall,  KJ 

U  specific  internal  energy  of  gas,  KJ/Kg 
R  gas  constant,  KJ/Kg  K 

Re  Reynolds  Number,  CmpD/v 

T  gas  temperature,  K3 

v  specific  volume,  rrw/Kg 

V  volume,  m^ 

X  dimensionless  pressure  defined  by  eqn.  (7) 
$  crankangle,  degree 

combustion  time  in  degree  crank 
e  engine  compression  ratio 

p  density,  Kg/m3 

y  dynamic  viscosity,  Kg/m  hr 

A  equivelance  ratio 

Subscripts 

a  at  the  beginning  of  compression  stroke 

f  fuel 

i  inlet  to  the  swirl  chamber 

m  main  chamber 

o  initial 

s  swirl  chamber 

t  tangential  port 

w  wall 


INTRODUCTION 

1 

The  design  of  internal  combustion  engines  is  necessarily 
becoming  more  scientific  as  the  standards  of  performance,  eco¬ 
nomy,  and  pollution  control  are  increased.  One  manifestation 
of  this  trend  is  the  increasing  use  of  the  mathematical  simu¬ 
lation  as  a  design  and  development  tool.  However,  an  engine 
design  based  on  a  mathematical  model  may  give  some  deviated 
results  due  to  the  assumptions  in  the  model.  Many  authors  jl, 
2,3  {.*  noticed  that,  along  with  other  information,  a  knowledge 
of  the  instantaneous  heat  transfer  to  the  cylinder  surface  is 
necessary  to  formulate  an  accurate  model  of  engine  processes. 

Many  authors  J  4—7  J  have  presented  expressions  for  the 
surface  heat  transfer  coefficient.  The  existing  formulae  are 
not  satisfactory  due  to  the  inaccuracy  in  predicting  the  gene¬ 
ral  magnitude  of  the  transfer  rate.  During  the  compression 
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stroke  heat  transfer  is  entirely  convective,  while  during  ex¬ 
pansion  stroke  it  is  convective  and  radiative.  Quite  large 
variations  in  general  magnitude  of  the  convective  component  of 
transfer  must  arise  from  engine  to  engine  or  even  from  point 
to  point  in  the  same  engine  because  of  variation  in  flow  pat¬ 
tern,  Theoretical  consideration  of  the  nonsteady  heat  trans¬ 
fer  was  explored  by  Pf laum  j 8  J  developed  by  Esler  J  9 {  and  ex¬ 
tended  by  Oguri  ] 10 J .  Assumptions  were  taken  in  these  analy¬ 
ses  for  the  turbulent  heat  transfer  near  the  surface.  In  ad¬ 
dition  correction  factor  was  applied  to  the  black  body  equa¬ 
tion  to  compensate  for  variation  from  black.  In  compression 
ignition  engines  combustion  proceeds  by  means  of  a  diffusion 
flame  created  in  unhomogeneous  fuel-air  mixture.  Thus,  the 
problem  involves  radiation  from  a  mixture  of  gases  at  high 
and  varying  temperature  to  complicated  wa3 1  shape.  These  facts 
clarified  that  it  is  impossible  to  calculate  accurately  the 
radiation  transfe  ;  during  the  combustion  expansion  period  and 
the  error  obtained,  when  using  correlations  based  on  experi¬ 
mental.  bases,  will  be  smaller. 

Numerous  investigations  have  been  carried  out  to  measure 
the  heat  transfer  rate  in  an  internal  combustion  engine. 

Nusselt  formula  Jllj  was  determined  on  the  basis  of  heat  loss 
measurements  from  spherical  bombs.  Brilling  ] 12 j  modified 
Nusselt' s  expression  by  introducing  a  term  for  the  speed. 
Eichelberg  J  4  J  and  Pflaum  {5]  conducted  experimental  work  to 
measure  the  instantaneous  heat  transfer  in  an  operating  engine. 
The  heat  flux  was  determined  by  measuring  the  wall  surface 
temperature.  They  used  it  as  a  boundary  condition  for  solving 
the  Fourier  unsteady  heat  transfer  equation  to  find  the  tem¬ 
perature  gradient  which  is  a  measure  of  the  heat  transfer  co¬ 
efficient.  Henien  j 13 |  reported  better  agreement  with  his 
measurements  when  replacing  mean  piston  speed  by  a  calculated 
instantaneous  gas  velocity.  Pischinger  ] 14 [  noticed  that  the 
radiant  transfer  cannot  be  neglected.  He  had  made  his  inves¬ 
tigation  on  several  diesel  engines  with  open  chamber  over  a 
speed  range  of  350-800  RPM.  Using  the  ratio  of  maximum  charge 
velocity  leaving  the  squish  area  to  piston  speed  he  was  able 
to  correlate  the  heat  fluxes.  It  was  found  that  the  correla¬ 
tion  varied  proportionally  to  the  ratio  of  the  charge  velocity 
to  the  piston  speed. 

Many  different  expressions  have  been  proposed  to  corre¬ 
late  the  surface  heat  fluxes  in  diesel  engines.  Attention 
was  focused  on  three  correlations  upon  which  most  other  corre¬ 
lations  have  been  based.  The  three  correlations  which  are 
considered  in  some  details  are  those  of  Nusselt,  Eichelberg, 
and  correlation  based  on  the  Reynolds  analogy  boundary  layer 
theory.  The  modifications  to  these  three  basic  forms,  which 
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have  been  suggested  by  various  authors  j  6 , 7  [are  also  con¬ 
sidered.  It  is  quite  clear  that  all  mentioned  formulae  fail 
in  one  way  or  another  to  meet  important  requirements.  Those 
which  are  dimensionally  nonhomogeneous  can  hardly  be  relied 
upon  for  extrapolation  to  conditions  different  from  those 
of  the  experiments.  Others  did  not  seen  to  allow  correctly 
either  for  unsteady  conditions  nor  for  swirl  chambers. 

Therefore,  study  of  the  instantaneous • surf ace. .heat 
transfer  in  an  indirect  injection  diesel  engine  has  been  car¬ 
ried.  An  experimental  plant  has  been  built  to  study  the  effect 
of  various  working  conditions  and  geometrical  parameters  of 
the  swirl  chamber  on  the  heat  transfer  coefficient  under  fir¬ 
ing  operation.  The  cylinder  head  was  prepared  with  a  special 
combustion  chamber  which  facilitated  the  variation  of  the  re¬ 
lative  swirl  chamhe1'  volume,  the  relative  connecting  channel 
area,  the  inclination  angle  of  the  connecting  channel,  and 
the  fuel  pressure.  Generalized  equation  is  established  for 
the  calculation  of  the  instantaneous  heat  transfer  to  swirl 
combustion  chamber. 


EXPERIMENTAL  INVESTIGATION 


In  view  of  the  requirement  for  a  detailed  understanding 
of  the  effect  of  various  factors  on  heat  transfer  in  swirl 
combustion  chambers,  an  experimental  set  up  has  been  built. 

The  eggine  was  a  4-stroke,  single  cylinder,  swirl  chamber 
diesel  engine  with  a  7C  mm  bore,  57  mm  stroke  and  a  compres¬ 
sion  fatio  17.  The  engine  was  coupled  to  an  electric  dyna¬ 
mometer  and  instrumented  for  obtaining  performance  data  such 
as  speed,  power  and  consumption  of  air  and  fuel.  The  general 
layout  of  the  experimental  plant  is  shown  in  Fig.  8.  The  swirl 
chamber  was  constructed  in  the  form  of  a  sphere  of  two  parts. 
The  two  parts  are  fixed  to  the  cylinder  head  by  two  stud  bolts 
on  the  top  corner  of  the  cylinder.  The  swirl  chamber  arrange¬ 
ment  is  shown  in  Big. 9.  The  lower  part  of  the  chamber  can  be 
changed  to  allow  for  the  variation  <5f  the  inclination  angle 
and  area  of  the  tangential  port.  The  swirl  chamber  can  be 
changed  also  for  the  variation  of  the  swirl  chamber  volume 
without  changing  the  compression  ratio.  The  swirl  chamber  was 
supplied  by  the  fuel  from  the  top  by  a  pintle  type  injector 
with  a  single  hole  of  diameter  1  mm  and  injection  angle  90 
degree . 

For  measuring  the  mean  surface  temperature  six  surface 
thermocouples  constructed  from  chrome l-Alumel  were  immersed 
deeply  near  to  the  surface  in  six  locations  in  the  periphery 
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FIG. 8  SCHEME  OF  EXPERIMENTAL  SET  UP  FOR  HEAT  TRANSFER 
STUDY 
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of  the  swirl  combustion  chamber.  The  fixation  of  thermo¬ 
couples  was  satisfied  such  that  the  disturbence  to  the  heat 
transfer  pattern  was  kept  to  a  minimum  and  the  thermocouple 
junction  temperature  was  considered  to  be  the  true  surface 
temperature.  A  selector  switch  and  a  potentiometer  were 
used  to  measure  the  time  averaged  value  of  the  inside  sur¬ 
face  temperature  of  the  chamber  wall. 

The  rapid  variation  of  gas  pressure  inside  the  com¬ 
bustion  chamber  together  with  the  crankangle  were  picked 
by  a  highly  sensitive  transducer  of  sernitivity  10. 95x10“ 
C/bar  and  natural  frequency  100  KHZ  and  an  external  X-gene- 
rator  consequently.  The  display  of  the  gas  pressure  was 
syncronized  with  the  crankshaft  by  a  disc  displaying  the 
signal  of  the  top  dead  centre.  A  strain  gauge  system  and  in¬ 
ductive  transducer  were  used  for  recording  fuel  line  pressure 
and  needle  lift  facilitating  the  calculation  of  the  instan¬ 
taneous  rate  of  fuel  injection.  A  sample  of  the  results  is 
shown  in  Fig  Jo.  A  precise  fuel  flow  meter  was  installed  bet¬ 
ween  the  tank  and  the  injector  to  measure  the  rate  of  fuel 
delivered  to  the  engine.  This  was  used  to  facilitate  the 
choice  of  the  suitable  discharge  coefficient  of  the  inject-  ^ 

The  data  recorded  during  an  engine  run  included  combus¬ 
tion  chamber  pressure,  crankangle,  top  dead  centre  and  six 
surface  temperatures.  Because  of  the  cycle  to  cycle  varia¬ 
tions  in  data  recorded  and  hence  in  surface  heat  flux,  it 
was  required  to  examine  a  large  quantity  of  data  if  meaning¬ 
ful  and  significant  results  were  to  be  obtained.  Thus,  a 
high  speed  multichannel  data  recording  and  processing  system 
was  developed.  The  output  from  each  transducer  was  magr ;f:'' 
and  recorded  by  an  8-channel  magnetic  tape  recorder.  A  ca.: 
bration  signal  was  recorded  at  the  start  of  recording.  There 
after  the  analog  signal  representing  chamber  pressure  was  dis 
played  back  from  the  tape  recorder  at  a  lower  speed  to  the 
analog-digital  converter  and  digitized  at  every  crank  for 
cycles.  The  output  signals  of  the  six  thermocouples,  fuel 
line  pressure  and  needle  lift  were  also  treated  in  the  same 
way.  The  output  of  the  50  cycles  was  then  processed  and  the 
average  cycle  pressure  variation  and  surface  temperature  were 
determined  and  recorded  onto  magnetic  tape  in  digital  form 
for  subsequent  processing  on  digital  computer  to  determine 
the  surface  heat  flux.  This  technique  facilitated  also  the 
attenuation  of  the  random  noise  introduced  by  the  electronic 
components . 


The  aim  of  the  experimental  investigation  was  the  study 
of  the  effect  of  volume,  throat  diameter,  throat  angle  of  the 
•wirl  combustion  chamber  on  heat  transfer  rate.  Thus,  the 


Fig. ( 10)  The  courses  of  swirl  chamber  gas  pressure,  fuel 

pressure,  needle  lift  &ihe  location  of  TDC  in 
relation  with  time . 
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essential  factors  were  varied  as  follows  : 

1-  Diameter  of  the  swirl  chamber  ranging  between  21.5  ram 
and  26  mm  satisfying  swirl  chamber  volume  ratios  ranging 
between  0.024  to  0.042. 

2-  Throat  diameter  of  tangential  port  ranging  between  7 
mm  to  10.5  mm  satisfying  throat  area  ratios  ranging  between 
0,1  to  0,15. 

3-  Throat  angle  of  swirl  chamber  ranging  between  45  and 
90  degree  satisfying  various  inlet  moment  of  momentum. 

4-  Engine  speed  ranging  between  1500  and  2500  RPM. 

5-  Injection  pressure  ranging  between  85  and  140  bar. 

6-  Engine  load  varying  between  1/4  and  full  load. 

One  engine  operating  condition,  which  is  given  in 
Table  1,  was  defined  as  the  standard  operating  condition. 
Any  variable  whose  value  is  not  specified  may  be  assumed  to 
be  the  standard  value. 


Table _ 1  Standard  operating  conditions 


Compression  ratio 

Engine  speed 

Fuel  injection  advance 

Inlet  temperature  and  pressure 

Inlet  valve  opens  and  closes 

Exhaust  valve  opens  and  closes 
Fuel 

Fuel  opening  pressure 
Throat  area  ratio 
Throat  angle 

Swirl  chamber  volume  ratio 
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2000  +  20  RPM* 

24+1  degree  CA  BTDC 
300  K,  1  bar 
704  degree  CA  and  216 
degree  CA. 

496  degree  CA  and  8deg.CA 
light  diesel  fuel. 

110  bar. 

0.12 

60  degree. 

0.024. 


HEAT  TRANSFER  MODEL  FOR  A  SWIRL  CHAMBER  DIESEL  gNGINE 

The  theory  used  here  simply  computes  the  heat  flux  from 
pressure  time  data.  The  mathematical  model  and  calculation 
procedure  is  thus  a  refined  model  of  obtaining  data  for  further 
correlation  studies. 

The  model  assumes  thermodynamic  equilibrium  at  each  instant. 
Thus  it  is  assumed  that  the  entire  cylinder  contains  homo¬ 
geneous  mixture  of  air  and  combustion  products  at  each  instant. 


Phenomena  such  as  fuel  vaporization,  mixing,  temperature 
gradients,  and  so  are  thus  ignored.  The  instantaneous  equi-  , 
librium  composition  are  considered  in  the  model.  With  these 
assumptions  the  energy  equation  during  the  compression  stroke, 
where  the  air  flows  from  the  main  combustion  chamber  to  the 
swirl  chamber,  becomes  ; 


d<{> 


dQ 

w 

d6 


dM.  6U  ..  6U_  dT  dM 

'  dTtUin+R  Tm)+  Ms  JT  3$  +  W  +Us“3? 


and  during  the  expansion  stroke  becomes 
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su 


6X 


s  dX 
d<j> 


6U  dT  dM 

s  s  s 

W  w  +Us“d$ 


(2) 

Before  fuel  burning  the  quantity  of  heat  released  QR  and  the 
term  (6U  /oX)  (dX/db)  are  equal  to  zero. 

The  internal  energy  of  the  products  of  combustion  of  fuel  and 
air  were  calculated  at  the  equilibrium  state  using  the  data 
from  JANAF  tables  }  14  j  .  The  fuel  is  assumed  to  take  the  form 
c16  h34  and  complete  combustion  is  assumed  due  to  presence  of 
excess  air  in  diesel  engines.  The  effect  of  dissociation  is 
largest  at  low  pressures,  high  temperatures  and  excess  air 
factor  near  unity.  However,  the  equilibrium  composition  was 
calculated  from  the  equation  of  combustion  for  various  tempe¬ 
ratures,  The  average  specific  heat  and  gas  constant  of  the 
constituents  were  calculated  and  were  then  fit  as  a  function 
of  T  P  to  achieve  a  great  save  of  computer  time.  The  final 
equation  for  the  internal  energy  and  the  gas  constant  were 
thus  found  to  be  given  by  ; 

U  =  (X-l)-a3  T  +  a4  T2  +  a5  T3  +  ag  (X-l)T 

for  X.  <1,6  and  T  >  1333  degree  K. 

U1=  U  +  14,2029 (T-1388.89)4'5371|l-10(X-l) }/P0*46024 


where , 

13) 

Const , 

al 

a2 

a3 

a4 

a6 

X  >  1.6 

T  <  1333K  - 

2888.78 

595.25 

0.668454 

2.713xlO-4 

-  7.40325xl0"S  0.25 

X  >  1.6 

T  >  1333K  - 

3466.2 

903.84 

0.65898  - 

3,8x10  4 

8.2329xlO~8  zero 

0.22751  +  0.063438X  -  3.33005xl0~3  X2  (4) 
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The  partial  derivatives  can  be  written  explicitly  in  terms 
of  temperature  and  equivelance  ratio.  The  time  derivative 
of  the  equation  of  state  gives  an  expression  for  dTs/d<{>  in 
the  following  form  : 


dTs/d*  =  Ts 


<3Ps/d$ 


dMs/dq> 


dR/d$ 


This  equation  is  used  to  eliminate  that  quantity  from  the 
energy  equations  1  and  2. 


The  inlet  mass  to  the  swirl  chamber  was  calculated  from 
the  model  suggested  by  ELKOTB  } 15 ] ,  in  which  a  relation  of 
the  instantaneous  pressures  inside  the  main  and  swirl  chambers 
is  derived  in  the  following  form  : 


d(v/vs) 

d« 


30  caAe’JVl 

"  *  vs 


3-n 

.  (X+v/v  -  \J  1-Xn_1 

where , 

X  =  (Pg/?)1/11 


(6) 

(7) 


P  =  j  Ce/as)  /  (X  +  v/vs) I"  (8)  * 


The  solution  of  this  differential  equation  gives  the 
instantaneous  /alue  of  the  dimensionless  swirl  chamber  pressure. 
The  swirl  chamber  pressure  values  are  calculated  from  a  smooth¬ 
ed  pressure-crankangle  using  a  second  degree  Lagrangian  in¬ 
terpolation  formula.  This  is  required  to  check  the  value  of 
the  discharge  coefficient  chosen  in  egn  6.  The  rate  of  mass 
flow  to  the  swirl  chamber  can  thus  be  written  in  the  following 
form  : 


This  is  equal  to  the  rate  of  increase  or  decrease  of  the  mass 
in  both  chambers. 

Before  combustion  proceeds  the  quantity  of  burned  fuel 
is  equal  to  zero  and  the  equivelance  ratio  in  both  chambers  is 
infinity.  After  ignition  the  quantity  of  fuel  increases  and 
the  equivelance  ratio  in  the  swirl  chamber  based  on  the  burn¬ 
ed  fuel  decreases.  An  increase  of  pressure  occurs  causing  a 
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flow  of  fuel-air  mixture  and  exhaust  gases  from  the  swirl 
chamber  to  the  main  chamber.  Therefore,  the  equivelance 
ratio  in  the  main  chamber  starts  to  decrease.  The  quantity 
of  burned  fraction  at  any  instance  from  the  beginning  of  com¬ 
bustion  was  estimated  from  the  equation  proposed  by  Setkijl6| 
as  follows  : 

•  mf(j/mf  =  1  "  lexP  “  6  *9  (4>/$c)m+1l 

where , 

m  =  (3.2  /  log  *c)  -1 

The  instantaneous  injected  fuel  quantity  was  determined  expe¬ 
rimentally  and  therefore  the  quantity  of  unburned  fuel  flow¬ 
ing  to  the  main  chamber  can  be  determined.  Thus,  the  quan¬ 
tity  of  heat  released  inside  the  swirl  chamber  is  calculated 
as  follows  : 

d  QR/d  4>  =  c.v  (d  mf/d<{>) 

For  small  interval  of  crankangle  the  equivelance  ratio 
at  the  end  of  the  .interval  is  calculated  from  the  calculated 
equivelance  ratio  at  the  beginning  of  the  interval  and  the 
rate  of  fueJ  burning  during  the  interval  d  <J>.  For  swirl  cham¬ 
ber  the  quantity  of  burned  fuel  at  the  end  of  interval  is  equal 
to  the  quantity  of  burned  fuel  at  the  beginning  plus 'the  burn¬ 
ed  fuel  during  the  interval  d<f>minus  the  quantity  of  burned 
fuel' flow  to  the  main  chamber.  Thus,  the  rate  of  change  of 
the  equiveiance  ratio  can  be  determined. 


EKSULTS  OF  THE  HEAT  TRANSFER  CALCULATIONS 

Calculations  of  instantaneous  heat  transfer  coefficient 
have  been  carried  out  using  pressure  crankangle  data.  The 
surface  temperature  of  the  swirl  chamber  must  be  estimated  i.rom 
experimental  data.  Therefore,  the  experimental  results  ob¬ 
tained  in  the  present  work  consists  mainly  of  measurement  of 
mean  wall  temperature,  measurement  of  pressure  time  diagram 
arid  measurement  of  fuel  delivery  law.  These  data  were  obtain¬ 
ed  from  a  swirl  chamber  air  cooled,  single  cylinder  research 
engine  at  various  operating  conditions  (engine  speed,  load  and 
fuel  opening  pressure)  and  various  geometrical  dimensions  of 
the  swirl  Chamber. 

The  effect  of  swirl  chamber  volume,  throat  area  of  the 
tangential  port,  inclination  angle  of  the  tangential  port  and 
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fuel  opening  pressures  on  the  mean  surface  temperature  was 
studied.  For  each  chamber  the  surface  temperatures  at  four 
locations  were  measured  at  various  speeds  while  the  factors 
were  kept  constants  at  the  standard  values  given  in  Table  1. 
The  mean  wall  temperature  was  calculated  from  the  arithmatic 
mean  value  of  the  four  measured  temperatures.  The  results  at 
various  engine  speeds  and  fuil  load  are  given  in  Figs.  (LI)  to 
(14),  One  notices  from  Fig.  (11)  that  the  mean  wall  temperature 
decreases  with  the  increase  of  engine  speed,  but  the  rate  of 
decrease  .is  relatively  great  in  swirl  chambers  with  greater 
diameters  than  with  small  diameters.  It  is  clear  also  that 
the  mean  wall  temperature  increases  when  the  diameter  of  the 
swirl  chamber  increases  as  the  heat  transfer  to  the  walls 
increases  with  the  increase  of  the  surface  area.  The  moment 
of  momentum  increases  also  with  the  increase  of  diameter  and 
consequently  the  gas  velocities  inside  the  chamber  increases 
giving  an  increase  to  the  heat  transfer  to  the  walls. 

One  notices  also  from  Fig.  ill}  that  the  mean  wall  tempe¬ 
rature  increases  when  the  throat  area  ratio  increases  or  dec- 
creases  than  0.1  since  the  cooled  air  at  the  entrance  of  the 
port  is  small  in  this  case  and  the  moment  of  momentum  in¬ 
creases  with  the  increase  of  the  throat  diameter.  The  varia¬ 
tion  of  the  inclination  angle  of  the  tangential  port  showed 
that  the  mean  wall  temperature  is  greater  in  swirl  chamber 
with  throat  inclination  of  45  and  90  degree  than  of  60  deg¬ 
ree,  as  shown  in  Fig.  (L3)  due  to  the  presence  of  cooled  air  at 
the  entrance  of  the  tangential  channel.  The  variation  of  the 
fuel  opening  pressure  shows  an  increase  of  the  mean  wall  tem¬ 
perature  with  the  increase  of  decrease  of  the  fuel  opening 
pressure  than  110  bar,  as  shown  in  Fig.  (14)  .  This  is  believed 
to  be  due  to  improvement  of  fuel  air  mixing  at  this  fuel  open¬ 
ing  pressure.  Similar  results  were  obtained  at  various  loads. 

Calculations  of  heat  release  in  these  swirl  chambers 
have  been  carried  out.  Samples  of  the  results  for  the  rate 
of  heat  release  are  shown  in  Fig.  (1^  .  The  first  peak  in  the 
heat  release  corresponds  to  the  period  of  rapid  pressure  rise 
and  lasts  for  only  about  2  degrees  C.A.  The  second  stage 
corresponds  to  the  period  of  gradually  decreasing  rate  of  heat 
release  and  lasts  about  13  degrees  C.A.  This  is  the  main  heat 
release  period.  Normally  about  80  percent  of  the  total  heat 
is  released  in  these  two  periods.  The  third  period  corres¬ 
ponds  to  the  tail  of  the  rate  of  heat  release  diagram  inwhich 
a  small  rate  of  heat  release  occurs. 

The  surface  heat  flux  and  the  heat  transfer  coefficient 
have  been  calculated  at  full  load  and  at  engine  speed  2000  RPM 
using  the  previous  calculated  heat  release,  mass  transfer  from 
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Fig  ii  delation  between  mean  wail  temperature 
and  engine  speed  for  different  swirl 
■  chamber  volume  ratios  at  full  load. 
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F3gC- 12  ) Relation  between  mean  wall  temperature 

and  engine  speed  for  different  throat  diamete* 
of  the  tangential  channel  at  full  load. 
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Fig(  13  delation  between  mean  wall  temperature 
and  engine  speed  for  differentinclination 
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Mean  wall  temperatures 
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Figc  14  delation  between  mean  wall  temperature 
and  engine  speed  for  different  opening 
pressures  of  the  injector  at  full  load. 


the  swirl  chamber  to  the  main  chamber.  Fig.  0-5)  ,  and  the  calcu¬ 
lated  gas  temperature.  The  instantaneous  gas  temperature  has 
been  calculated  using  the  equation  of  state  for  which  the  in¬ 
stantaneous  gas  constituents  with  dissociation  were  calculated. 
The  calculated  heat  flux  and  the  heat  transfer  coefficient  for 
various  swirl  chamber  diameters  are  shown  in  Fig.  (L6)  .  The 
surface  heat  flux  increases  firstly  reaching  maximum  value  at 
the  T.D.C.  and  decreases  again.  The  values  of  the  surface 
heat  flux  increases  with  the  increase  of  the  swirl  chamber 
diameter.  The  moment  of  maximum  heat  flux  is  not  affected  with 
the  change  of  chamber  volume.  The  heat  transfer  coefficient 
increases  also  similarly,  as  the  surface  heat  flux,  with  the  in 
crease  of  chamber  volume  and  the  maximum  value  of  the  heat  tran 
sfer  coefficient  is  nearly  at  the  T.D.C. 

The  effect  of  the  throat  diameter  and  the  inclination 
angle  of  the  tangential  port  on  the  surface  heat  flux  and  the 
heat  transfer  coefficient  have  been  studied  and  the  results 
are  illustrated  in  Figs.  (17)  and  (18).  The  influence  of  the 
throat  diameter  and  the  inclination  angle  is  not  great  compar¬ 
ed  to  the  effect  of  the  chamber  volume.  Minimum  value  of  heat 
transfer  coefficient  is  obtained  when  the  throat  diameter  was 
8.5  nun  and  the  inclination  angle  60  degrees.  The  location  of 
the  maximum  value  is  not  affected  with  the  change  of  the  throat 
diameter  and  the  inclination  angle.  The  effect  of  fuel  open¬ 
ing  pressure  has  been  also  studied  and  the  results  are  illus¬ 
trated  in  Fig.  (19).  it  was  found  that  the  heat  transfer  co¬ 
efficient  has  greater  values  with  fuel  opening  pressure  of  85 
bar  .and  smaller  values  with  fuel  opening  pressure  of  110  bar. 
This'  means  that  when  the  fuel  opening  pressure  increases  or 
decreases  than  110  bar  the  quantity  of  fuel  burned  in  the  swirl 
chamber  increases. 


CORRELATION  AND  COMPARISON  OF  HEAT  TRANSFER  COEFFICIENT 

It  is  clear  from  the  peevious  discussion  that  the  heat 
transfer  process  is  a  very  complex  phenomesan  and  the  early 
empirical  realations  were  uns?»tisfactory  due  to  their  suitabi¬ 
lity  to  specified  engine.  An  obvious  solution  to  this  compli¬ 
cated  problem  is  to  derive  the  dimensionless  groups  assuming 
steady  state  system  at  any  instant.  The  variables  affecting 
the  heat  transfer  process  in  this  case  are  as  follows  ; 

h  =  f(D,  ds,  dt,  C^,  Cp,  yi,  K,  PfQ) 

The  application  of  dimensionless  analysis  reduces  the  variables 
to  the  following  dimensionless  groups  : 


P,  *  110  bar 


rfj  »  8-5mm. 

d  *  60 


COMPARISON  BETWEEN  EXPERIMENTAL  WORK  AND  EM 
PIRICAL  CORRELATION 
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and  the  equation  of  the  heat  transfer  becomes  : 

Nu  =  A  ReaPrb  EuC  (dg/dfc) e 

Using  the  experimental  results  the  indices  of  the  dimension¬ 
less  quantities  were  determined  by  plotting  the  Nusselt  num¬ 
ber  versus  one  of  these  quantities.  The  indices  that  fit  the 
experimental  values  by  a  smooth  curve  are  the  required  values 
for  the  generalized  equation.  Thereafter,  plotting  the  Nusselt 
number  versus  the  dimensionless  quantities  raised  to  .the  ob¬ 
tained  indices  as  shown  in  Fig.  (2o) »  the  constant  of  the  gene¬ 
ralized  relation  was  obtained.  The  generalized  relation  has 
been  found  in  the  following  form  : 


Nu  ®  1.175xlO~4  Re0*5  P 


PA  Eu-1  ’ 5  (dg/ dj.) 1  * 2 


Comparing  the  results  of  this  work  with  other  previous 
studies  such  as  Woschni,  Annad  (1970),  Annad  (1963)  and  it  is 
clear  from  Fig. (21)  fchat  Nusselt  number  in  this  study  for  the 
same  Reynolds  number  is  greater  than  all  other  values.  This 
was  expected  as  the  velocity  and  turbulence  of  the  mixture  in 
the  swirl  chamber  is  more  than  in  the  open  chamber. 
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CHAPTER  5 


PREDICTION  AND  MEASUREMENT  OF  FLOW  AND  HEAT 
TRANSFER  IN  MOTORED  DIESEL  ENGINES  SWIRL  CHAMBERS 


SUMMARY 


This  paper  is  concerned  with  the  prediction  and  measure¬ 
ment  of  flow  and  heat  transfer  in  motored  diesel  engine  swirl 
chambers . 

Measurements  are  mainly  carried  out  by  a  temperature- 
compensated  hot-wire  anemometer  and  a  high  speed  camera.  Pre- 
ditions  are  based  on  the  solution  of  the  finite-difference 
form  of  the  governing  differential  equations  for  the  transport 
of  mass,  momentum  and  energy  by  way  of  digital  computers.  In 
general,  the  obtained  agreement  between  predictions  and  meas¬ 
urements  is  fairly  good  which  confirms  the  accuracy  of  the  des¬ 
cribed  prediction  procedure. 


NOMENCLATURE 

Ac  influence  coefficient 

Ap  equal  to  £AC 

Cp,Ci,C2  Constant  of  the  turbulence  model 

Cpm  mean  piston  velocity 

h  total  enthalpy 

K  turbulence  kinetic  energy 

mc  mass  flow  rates  across  cell  boundaries 

Mft  equals  to  (pV)g/ot  old  values  at  time 

Mn  equals  to  (pV)pn/6t  new  values  at  time  t+6t 

Pp  pressure 

p‘  pressure  correction 

u  tangential  velocity 

u  velocity  vector  at  any  point 

r  radial  position 

Sp  coefficient  of  linearized  source  term 

Su  coefficient  of  linearized  source  term 

S  source/sink  term  for  variable  4> 
t  time 

v  *  radial  velocity 

V  cell  volume 

0  equal  to  (3  p  /3p)T 

If  effective  diffusivity  coefficient  for  variable  $ 
e  dissipative  rate  of  turbulent  kinetic  energy 

0  angular  position 

$  arbitrary  dependent  variable 
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y  laminar  viscosity 

yeff  effective  viscosity 

p  density 

ahfak'aF  constants  of  turbulence  model  for  h,  K  &  e. 

£c  summation  for  neighbouring  nodes  of  typical  grid  node. 

Subscripts. 

E,W,N,S  East,  West,  North  and  South  neighbouring  nodes 
e,w,n,s  midway  nodes  between  each  node. 


INTRODUCTION 

The  availability  of  very  large,  high  speed  digital  com¬ 
puter  has  encouraged  efforts  to  simulate  the  compression  igni¬ 
tion  engines.  Complete  mathematical  model  for  the  calcula¬ 
tion  of  the  heat  release  in  a  swirl  chamber  of  a  diesel  en¬ 
gine  is  still  not  available.  The  heat  release  in  multifuel 
engines  depends  to  a  great  extent  on  the  hydrodynamic  mixing 
and  transport  processes,  the  fuel  atomization  and  evaporation, 
the  chemical  reaction  mechanism  and  reaction  rates,  heat  tran¬ 
sfer  and  environmental  conditions,  all  of  which  vary  inside 
the  swirl  chamber.  The  air-fuel  mixing  process  perior  to 
autoignition  is  important  in  diesel  engine,  in  general  and  in 
swirl  chamber  diesel  engine  in  particular.  Moreover  the  mix¬ 
ing  process  controls  combustion  which  in  turn  controls  the 
rate  of  pressure  rise,  noise  emission  and  pollutant 
formation  in  engine  cylinders  j 1 , 2 , 3  J . *  The  study  of  pure 
phenomena  in  subclasses  will  help  the  solution  of  the  more 
complex  combustion  phenomena  resulting  from  multifuel  combus¬ 
tion. 

Therefore,  a  mathematical  prediction  model  capable  of 
predicting  the  air  velocity  inside  swirl  combustion  chamber 
is  required  for  the  development  of  multifuel  engines  as  well 
as  saving  time  and  cost. 

There  have  been  several  experimental  investigations,  as 
a  consequence  of  the  development  of  swirl  chambers;  for  the 
determination  of  the  swirl  chamber  parameters.  However,  it 
was  found  that  the  study  of  flow  is  the  more  accurate  method 
for  determination  of  swirl  chamber  parameters.  Theoretical 
investigations  have  been  therefore  reported  | 4 , 5 , 6  f  .  Most 
of  these  investigations  are  based  on  the  evaluation  of  the 
intensity  of  charge  motion  which  relies  in  emperical  correla¬ 
tions,  Moreover,  they  are  based  on  the  application  of  evaluat¬ 
ing  parameters  recommended  from  the  analysis  of  the  present 

*  A  list  of  references  is  given  at  the  end  of 
the  chapter . 


swirl  chambers.  Such  evaluating  parameters  have  a  quite 
wide  range  that  they  cannot  give  sufficient  accurate  results. 
Recently,  ELKOTB  J7j  proposed  a  new  theoretical  model  based 
on  the  equality  of  the  moment  of  momentum  resulting  from  the 
ejection  velocity  to  the  swirl  chamber  and  the  moment  of  mo¬ 
mentum  of  the  rotating  air  inside  the  swirl  chamber  provid¬ 
ed  that  the  air  motion  in  the  central  zone  is  solid  vortex 
and  that  in  periphery  is  free  vortex.  This  model  can  define 
the  average  velocity  to  the  swirl  chamber  and  the  moment  of 
momentum  of  the  rotating  air  inside  the  swirl  chamber  pro¬ 
vided  that  the -air  motion  in  the  central  zone  is  solid  vor¬ 
tex  and  that  in  Periphery  is  free  vortex.  This  model  can 
define  the  average  velocity  and  can  not  facilitate  the  study 
of  the  interaction  between  fuel  and  air.  Over  the  past  de¬ 
cade  a  fairly  considerable  effort  has  been  expended  on  com¬ 
putational  procedures  with  the  aid  of  mathematical  models, 
for  reciprocating  engines  j  8 , 9 , 10 { .  These  prediction  proce¬ 
dures  are  based  mainly  on  the  solution  of  the  finite  dif¬ 
ference  for  the  governing  differential  equations  for  the 
transport  of  mass  momentum  and  energy. 

Although  these  attempts  have  been  made  to  ipredict  the 
properties  of  turbulent  motion  by  mathematical  models,  no 
explicit  consideration  of  the  physical  properties  in  swirl 
chamber  is  taken.  These  includes  the  continuous  variation  of 
the  medium  pressure  and  temperature  and  heat  transfer  during 
flow.  Moreover,  the  shape  of  the  combustion  chamber  with  its 
effect  on  the  air  flow  was  not  included. 

,  The  objective  of  this  work  is  to  develop  a  complete 
mathematical  model  of  the  prediction  of  the  turbulent  flow 
field  in  swirl  chambers.  Experimental  measurement  of  the  air 
flow  velocity  was  planned  to  satisfy  careful  comparison  over 
a  wide  range  with  the  predicted  values  to  obtain  reliable  in¬ 
formation  about  the  constants  of  turbulence  model.  Actually 
these  constants  depend  on  local  turbulence  values.  This 
proposed  prediction  procedure  is  required  to  investigate 
speedily  and  with  reasonable  detail  and  accuracy  the  combus¬ 
tion  of  fuel  with  various  chemical  properties  as  well  as  the 
effect  of  various  factors  on  the  combustion  process  needed 
for  the  development  of  multifuel  engines. 


EXPERIMENTAL  INVESTIGATION 

A  careful  comparison  over  a  wide  range  of  data  can  give, 
in  addition  to  the  mean  velocity  field,  reliable  information 
about  the  constants  of  turbulence  model.  Therefore  an  experi¬ 
mental  setup  has  been  built  up  to  facilitate  measurements  of 
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the  air  velocity  components  inside  the  swirl  chamber  with  the 
highest  possible  accuracy.  The  experimental  swirl  chamber 
has  been  built  up  in  one  of  the  engine  cylinders  of  a  2-cy¬ 
linder,  4-stroke,  water  cooled  medium  speed  diesel  engine. 

The  other  cylinder  has  been  used  for  motoring  the  engine  to 
the  required  regime.  The  engine  was  directly  coupled  with 
hydraulic  dynamometer  for  engine  loading  and  equipped  with 
various  instruments  required  for  measuring  the  main  engine 
parameters. 

The  scope  of  this  investigation  includes  the  measuring 
of  the  instantaneous  velocity  components  and  finding  out  the 
effect  of  various  constructional  and  working  conditions  on 
the  air  flow  inside  swirl  chambers.  For  this  reason  the  cy¬ 
linder  head  was  modified  and  a  swirl  chamber  was  constructed 
and  attached  instead  of  the  original  chamber.  A  Schematic 
drawing  of  the  experimental  setup  is  shown  in  Fig. (22).  A 
cylinderical  combustion  chamber  of  48.5  mm  in  diameter,  Fig. 

(23)  satisfying  a  relative  swirl  chamber  volume  ratio  of  0.507 
was  constructed  with  an  optically  flat  fused  quartz  window 
on  one  of  its  sides.  The  swirl  chamber  was  built  up  with  com¬ 
pression  ratio  17.  Two  oblique  washers  were  fitted  on  either 
sides  of  the  quartz  plate  to  avoid  stress  concentration.  The 
quartz  lens  was  fastened  by  a  self  locking  screwed  collar 
facilitating  easy  firm  and  quick  assembly.  The  test  swirl 
chamber  was  fitted  with  a  cooling  jacket  supplied  with  gly- 
srine  from  an  independent  circulating  system  to  control  its 
surface  temperatnre.  Connecting  ports  with  various  confi¬ 
gurations  were  constructed  to  supply  air  to  the  swirl  chamber 
from  the  main  chamber.  The  chamber  was  designed  such  that 
variation  of  the  relative  swirl  volume  ratio [7 |,  the  connect¬ 
ing  port  area  ratio  J  7  J  and  the  connecting  port  inclination 
angle  are  possible.  The  intake  and  exhaust  manifolds  of  the 
test  cylinder  were  separated  from  those  of  the  second  motor¬ 
ing  cylinder. 

* 

The  chamber  was  equipped  with  a  temperature  compen-  » 

sated  hot  wire  anemometer  to  measure  the  air  velocity  compon¬ 
ents,  The  probe  of  the  hot  wire  was  connected  to  an  ultra- 
voilet  oscillograph  through  the  compensating  bridge.  It  was 
equipped  also  with  two  f errous-constantan  thermocouples  for 
measuring  the  surface  temperature.  The  quartz  window  was 
used  to  photograph  the  air  motion  inside  the  swirl  chamber 
with  a  high  speed  camera.  The  gas  pressure,  variations  in¬ 
side  the  main  and  swirl  combustion  chambers  together  with  the 
crankangle  were  picked  by  a  highly  sensitive  Piezo-electric 
transducers  and  an  external  X-generator.  The  display  of  the 
gas  pressure  was  syncronised  with  the  disc  displaying  the  sig¬ 
nal  of  the  top  dead  centre.  The  data  recorded  during  an  engine 


Fig (  22  }  Scheme  Of  E: 
Measurement 


FIG. 23  SCHEME  OF  THE  TEST  COMBUSTION  CHAMBER  FOR 
AIR  FLOW  MEASUREMENTS. 
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run  included  combustion  chamber  pressure,  velocity  variation 
for  two  inclination  angles  of  the  hot  wire  at  each  location, 
crank-angle,  top  dead  centre  and  two  surface  thermocouples. 
Because  of  the  cyle  to  cycle  variations  in  data  recorded  and 
hence  in  air  velocity  components,  it  was  required  to  examine 
a  large  quantity  of  data  to  obtain  a  quite  sufficient  re¬ 
asonable  results.  Thus,  a  high  speed  multi  channel  data  re¬ 
cording  and  processing  system  was  applied.  The  output  from 
the  hot  wire  anemometer  was  magnified  and  recorded  by  a  mag¬ 
netic  taperecorder .  A  calibration  signal  was  recorded  at 
the  start  of  recording.  Thereafter,  the  analog  signal  re¬ 
presenting  the  air  velocity  variation  at  two  inclination 
angles  were  played  back  at  a  lower  speed  to  the  analog/digi¬ 
tal  converter  and  digitzied  at  every  crank  angle  for  50 
cycles.  The  output  of  the  50  cycles  was  then  processed  and 
an  average  cycle  of  the  velocity  variations  was  determined 
and  recorded  on  magnetic  tape  in  digital  form  for  subsequent 
processing  on  digital  computer  to  determine  the  mean  velocity 
and  the  random  velocity. 


THE  GOVERNING  CONSERVATION  EQUATIONS 


The  flow  inside  diesel  engines  swirl  chambers  is  a 
transient  turbulent  flow.  For  the  case  of  motored  engines, 
thermal  radiation  may  be  neglected.  The  flow  is,  thus,  go¬ 
verned  by  the  differential  conservation  equations  of  energy, 
mass  and  momentum.  In  the  present  investigation,  the  swirl 
chamber  is  cylindrical  and  the  flow  may  thus  be  assumed  two- 
dimensional. 

The  equations  governing  the  flow  in  a  cylindrical  swirl 
chamber  of  a  diesel  engine  can  be  obtained  by  decomposing  the 
velocity,  pressure,  density  and  enthalpy  into  a  mean  and  a 
fluctuating  value  respectively  [ 9 , 11 [ . 

The  turbulent  diffusion  fluxes  which  appear  in  the  go¬ 
verning  equations  are  modeled  using  Boussinesq  approach  in 
terms  of  the  mean-flow  gradients  and  eddy  dif fusivities . 

These  eddy  dif fusivities  are  determined  by  solving  addition¬ 
al  differential  conservation  equations  for  the  time-averaged 
kinetic  energy  of  turbulence  k  and  its  dissipation  rate  e. 
This  kind  of  model  was  originally  developed  for  steady  tur¬ 
bulent  flows  }12|  and  has  been  modified,  to  introduce  com¬ 
pressibility  effects  by  analogy  with  the  laminar  stress  ten¬ 
sor  ,  by  RAMOS  J  9  J  . 
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Within  the  above  framework,  the  governing  set  of  equa¬ 
tions  may  be  compactly  represented  in  terms  of  a  single  gene¬ 
ral  equation  for  an  arbitrary  dependent  variable  4>  ; 

It  (p*)  +  r  I  {rpv*)+  k  Ig  (pu*>-  F  h  irTM] 

~  hh  (1i  ft*  =  s<}»  <a> 

where  $  stands  for  radial  velocity  v,  tangential  velocity  u, 
total  enthalpy  h,  k  and  e,  t  is  time,  p  is  density,  r^and 
are  respectively  the  effective  diffusivity  coefficients  and 
source/sink  terms  for  variable  <f>  .  Equation  (1)  represents 
also  the  continuity  equation  by  replacing  <J>  by  1.  The  defi¬ 
nition  of  r<j,  and  are  given  in  table  (1).  In  table  (l)yef 
is  the  effective  viscosity  which  is  given  by  ; 

yeff  =  p  +  Cpp  k2/£  C2) 

where, 

y  is  the  laminar  vis*  .ty,  ,  Cj_,  C2,  crh,  and 

ae,  are  constants  of  the  turbulence  model  and  are  given  in 
ref,  J  9 , 11 1  $  is  the  velocity  vector  at  any  point  in  the 

swirl  chamber. 

Table  (1):  Definitions  of  I’^and 
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4> 

K 


yeff/ak  G  "  p£ 


yeff^ae  K  tciG_C2pe^ 


_  i0/l  3u  ,  v.2  ~,3v,2  ,1  3 v  3u 

G~  yeffl2<r  3  0  +  r  2(Sr)  (r  3  0  3r 


|)2!-  |  div(u) 


(heff  div(uj  +  pK) 


BOUNDARY  CONDITIONS 


The  imposition  of  the  boundary  conditions  for  diesel  en¬ 
gine  swirl  chambers  is  not  as  straight-forward  as  it  may  seem. 
At  chamber  walls  :  the  velocity  components  all  obey  the  no 
slip  conditions;  the  walls  are  insulated  and  hence  the  wall 
heat  flux  is  taken  equal  to  zero  in  this  stage;  the  turbulence 
fluctuations  and  their  dissipation  rate  are  also  zero.  Wall 
functions  are  used  for  grid  nodes  adjacent  to  the  chamber  walls 
to  avoid  using  many  grid  lines  in  the  boundary  sub-layer  1 13  j . 
At  the  inlet  plane,  i.e.  the  tangential  port  which  connects 
the  swirl  chamber  with  the  main  chamber,  the  inlet  velocity, 
density  and  temperature  vary  with  crank  angle  Cor  time) .  They 
are  specified  by  solving  an  ordinary  differential  equation  for 
the  flow  from  the  main  chamber  to  the  swirl  chamber,  or  vise 
versa,  through  the  tangential  port,  in  the  manner  explained 
in  ref.  J 14  J  . 


THE  FINITE-DIFFERENCE  EQUATIONS 


For  the  purpose  of  deriving  the  finite-difference  equa¬ 
tions,  the  swirl  chamber  is  overlaid  with  a  grid  of  nodes,  for¬ 
med  by  the  intersections  of  meridunal  lines  and  circles  (Fig. (24) 
at  which  uhe  scalar  variables  h,K,e  and  P  are  stored;  while 
the  velocities  are  located  mid-way  betv/een  the  pressures  wh^ch 
drive  them.  Equation  (1)  is  formally  integrated  over  imagin¬ 
ary  control  volume  or  'cell'  surrounding  each  variable  and  for 
a  time  increment  5t.  The  finite-difference  approximation  to 
equation  (1)  is  written  as  ; 


VMp'Spl4,p'’  =  Ac  C  +  Mp  +p  +  su  (3> 


where , 

A  =  Ic  Aq,Ic  denotes  summation  over  the  four  neigh- 
bouringpnodes  or  a  typical  grid  node  P,c  denotes  neighbouring 
node  of  a  grid  node  P,  ’o’  and  ’n'  denote  'old*  and  'new'  and 
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'new'  values  at  times  t  and  t+ot  respectively,  Mp  equals (pV)p/ 
6t,  V  is  cell  volume,  Sg  and  Su  are  coefficients  of  a  lineari¬ 
zed  source  term  obtained  by  integrating  S^.hc  are  influence 
coefficients  which  give  the  combined  effect  of  diffusion  and 
convection. 


The  continuity  equation  is  integrated  to  give  ; 


Mn  -  M°  +  I  m 
p  p  c  c 


where , . 


0 


(4) 


mc  are  mass  flow  rates  across  cell  boundaries.  Equation 
(3)  can  be  solved  for  u,v,h,k  and  e  in  space  and  time  domains 
if  the  pressure  field  is  known  at  each  time  step.  The  solu¬ 
tion  of  Equation  (3)  may  be  obtained  economically  if  coeffi¬ 
cients  Aj., Sp  and  Su  are  computed  at  'Old'  time  t  and  leaving 
sjip  and  to  be  calculated  in  the  'new*  time  t+6t.  This  so¬ 
lution  approach  parallels  that  of  a  parabolic  flow  which  is 
justifiable  because  equation  (1)  is  indeed  parabolic  in  time 
Jlbj.  The  set  of  difference  equation  (3)  are  linked  through 
pressure  differences  terms. 


OUTLINE  OF  THE  SOLUTION  PROCEDURE 


At  any  instant  of  time  t  the  ©°  fields  of  the  variables 
are  known,  predictions  for  a  time  increment  ot  are  then  ob¬ 
tained  by  solving  the  difference  equations  for  the  ©n  fields 
by  a  marching-integration  algorithm.  The  latter  then  becomes 
the  <ti0  for  the  prediction  of  the  next  time  step;  in  this  way 
the  solution  is  marched  forward  in  time  until  the  desired  per¬ 
iod  has  been  covered.  The  next  section  is  devoted  for  the  so¬ 
lution  of  the  difference  equations. 


SOLUTION  OF  DIFFERENCE  EQUATIONS 

The  energy  equation  is  first  solved,  using  the  old  fields 
of  h,u,v,p  and  (3P/3t),  to  yield  the  new  enthalpy  field  (hn} 
and  hence  the  new  temperature  field  (Tn) . 

The  next  step  is  to  obtain  a  preliminary  set  of  new  ve¬ 
locity  components  v*n  and  u*n  for  an  estimated  new  pressure 
field  P  n  by  solving  their  difference  equations;  P*n  is  actual¬ 
ly  estimated  as  equal  to  the  'old*  pressure  P°. 

The  computed  velocities  v*n  and  u*n  will  not,  in  general, 
>  satisfy  the  local  continuity  equation  (4)  but  will  produce  a 
net  mass  source  at  each  grid  node.  The  pressure,  density  and 
velocities  are  then  corrected  so  as  to  reduce  the  mass  source 
at  each  grid  node  to  zero  in  the  following  manner  ; 
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where , 

the  *  indicates  guessed  or  preliminary  value,  P*  is  pres¬ 
sure  correction,  3  is  defined  as  (8p/9P)T,  Dv  and  Du  are  eva¬ 
luated  from  the  relevant  momentum  equation,  e.g.Du  =  u"/ (PD-Pg) , 
subscripts  N  and  E  denote  neighbouring  nodes,  as  shown' in 
Fig.  (24)*. 


Equations  (6-8)  are  substituted  into  the  continuity 
equation  (4)  to  give  ; 

(Ap-Sp)  Pp  =  Zc  Ac  Pc  +  Su  {9) 

where,  *  * 

A  =  £  A  ,  S.,=-  (Md  -  Mp  -5-  rcmc)  is  the  local  continuity 

imbalance  based  on  u  *  v  and  p  ;  S0  and  Ap  are  coefficients 

which  are  defined  '  as  follows  ; 


★  *n 

Sp  =  “VpV5t  '  £c(8cmc/pc  } 


(10). 


*  n  4  x  ■***  **  * 

AE  =  pe  aeDue  “°'5  8eVpe  ;Aw=  pw  W0'5  8VpW 


t  * 


=  °n  anDvn-°-5*nmn/pn  ;  AS  *  Ds  asDvs+0-5Ssras/p£ 


where , 


subscripts  E,T,v,N,S  and  e,  w,n,s  a.re  neighbouring  nodes 


and  corresponding  nodes  mid-way  between  each  node  and  a  typi¬ 
cal  nodes  P  respectively  (Fig .(24).  In  proving  equation  (9)  terms 
involve  P ’ ^  were  ignored  because  P*  itself  should  be  small. 


It  should  be  noted  that  the  pressure  correction  equation  dees 
not  include  P*  of  ‘old’  time  because  there  is  no  question  of 


correcting  'old'  velocities  and  pressure  fields. 


Equation  (9)  is  solved  to  yeild  the  P*  field  which  is 
required  to  correct  P,6,v  and  u  fields.  The  new  K  and  e  are 
then  computed  from  their  difference  equation  (3) . 
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OTHER  FEATURES 


The  difference  equations  (3)  &  (.9)  are  modified  at  cells 
adjoining  the  swirl  chamber  boundaries  to  incorporate  the  con¬ 
ditions  imposed  there.  Ir.  addition  to  P ' ,  a  global  pressure 
correction  P'  is  computed  based  on  overall  continuity  balance 
and  is  added  to  the  existing  pressure  field  before  solving  the 
pressure  correction  equation  (11) .  The  method  of  solution 
used  to  solve  equations  (3)  &  (9)  is  a  Gauss-elimination  line- 
by-line  double-sweep  technique.  While  equation  (3)  converges 
very  fast,  equation  (9)  needs  mere  sweeps  and  some  times  more 
than  one  iteration  to  reduce  the  residual  mass  errors  to  an 
acceptable  level.  The  modified  form  of  the  pressure  correc¬ 
tion  equation,  developed  here,  gives  a  fast  convergence  which 
reduces  the  total  computational  time.  The  time  increment  fit 
which  is  used  in  the  present  work  is  2.22x10”^  s  which  is 
equivelent  to  2°  crank  angle. 


PREDICTION  AND  MEASUREMENT  RESULTS 

The  com-'>utational  grid  inside  the  swirl  chamber  consists 
of  10  angular  and  8  radial  locations.  Finer  grid  could  be 
employed,  at  _he  expense  of  an  increase  in  computational  time; 
however  the  increase  of  number  of  grid  nodes  is  not  warranted 
at  this  stage.  The  computation  is  started  at  zero  time,  de¬ 
fined  as  the  time  at  which  the  compression  stroke  has ‘started. 
The  inlet  values  to  the  swirl  chamber  are  taken  from  the  pre¬ 
vious'  computational  method,  reported  in  } 15 [  .  Computations 
are  carried  out  for  a  diesel  engine  with  a  swirl  chamber  of 
48.5  mm  diameter,  tangential  port  of  diameter  15  mm  and  inclin¬ 
ed  angle  45°  to  the  cylinder  axis,  engine  speed  1500  RPM, 
and  compression  ratio  17.  The  computations  are  carried  at  2° 
crank  angle  intervals.  The  resulting  velocity  field  computed 
at  various  times  is  shown  in  Figs  (25)and  (26)a.n<?  the  corresponding 
contours  of  the  turbulence  intensity  distributions  (  V?K/3"/  Cpm, 
where  Cpin  is  the  mean  piston  velocity  (C  =  7m/s)  is  shown  in 
Fig. (27). 

In  time  sequence,  one  can  observe  the  features  of  the 
flow  evolution.  At  crank  angle  18°,  it  can  be  seen  that  the 
air  jet  exists  at  the  tangential  port  with  a  small  velocity 
which  is  principally  directed  from  both  sides  of  the  port  to 
the  opposite  wall  surface.  The  flow  in  the  swirl  chamber  dur¬ 
ing  the  initial  part  of  the  compression  stroke  is  almost  pure 
radial  flow.  After  some  time  for  example  at  58°  crank  angle, 
the  pressure  gradient  that  have  built  up  at  the  opposite  sur¬ 
face  have  already  initiated  rotational  flow  in  the  bulk  of  the 
swirl  chamber.  Beforu  the  initiation  of  the  complete  rotation¬ 
al  flow  a  vortex  is  noticed  at  the  right-hand  side  of  the  jet. 
Before  the  end  of  the  compression  stroke  the  radial  pressure 
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gradients  have  strengthened  to  the  point  of  inducing  strong 
flow  recirculation  which  was  observed  from  the  high  speed  pho¬ 
tographing  of  the  flow  field  as  shown  in  Fig.  (28)  .The  strong 
recirculating  flow  inside  the  swirl  chamber  produces  a  large 
tangential  velocity  component  along  the  chamber  surface  and 
small  tangential  velocity  component  at  the  centre,  as  shown 
in  Fig.  (29). It  is  noticed  also  that  the  tangential  velocity 
component  increases  during  compression  stroke  reaching  its 
maximum  value  at  about  15°  crank  BTDC,  as  shown  in  Fig.  (30). 

At  the  begining  of  the  expansion  stroke  the  recirculating  flow 
is  unable  to  reverse  its  direction  to  flow  through  the  tangen¬ 
tial  port  in  spite  of  the  absence  of  a  jet  flow  into  the  swirl 
chamber,  as  shown  in  Fig.  (26)  .Before  the  complete  reverse  of 
the  flow  a  vortex  is  noticed  at  the  opposite  side  to  the  port. 

It  should  be  mentioned  that  by  this  time,  the  flow  velocities 
are  low. 

The  turbulence  contours  during  compression  stroke  are 
also  shown  in  Fig  .(27)  in  which  it  is  noticed  that  the  turbu¬ 
lence  levels  at  the  begining  of  the  compression  stroke  are  low. 
By  increase  of  time  the  inflow  to  the  swirl  chamber  increases 
and  separates  at  the  begining  near  to  the  tangential  port  to 
form  eddies  at  either  sides  of  the  jet.  Turbulence  int^sity 
generation  is  particularly  strong  reaching  maximum  value  of 
1.33  near  to  the  entry  port.  Substantial  decrease  of  the  tur¬ 
bulence  intensity  has  been  obtained  by  BTDC  where  a  complete 
recirculating  flow  is  established.  This  decrease  of. turbu¬ 
lence  intensity  presumably  being  a  consequence  of  the  gradual 
reduction  of  the  inflow  rate  by  12°  crank  angle  as  the  piston 
apprbaches  the  end  of  the  compression  stroke. 

Figure  (31 )  shows  the  variation  of  the  turbulence  in¬ 
tensity  along  the  radius  (plane  2) .  Maximum  turbulence  inten¬ 
sity  is  noticed  near  chamber  walls  while  the  minimum  intensity 
occurs  at  a  radius  of  about  1.0  mm.  Turbulence  irtensity  va¬ 
riation,  at  10  mm.  from  chamber  walls  with  crank  angle  is  depec- 
ted  in  Fig. (32).  Turbulence  intensity  is  very  low  at  the  begin¬ 
ing  of  compression.  It  then  increases  rapidly  during  the  com¬ 
pression  stroke  reaching  a  maximum  valae  at  the  end  of  com¬ 
pression  where  the  turbulence  intensity  starts  to  decrease 
again  during  the  expansion  stroke. 

Comparisons  of  measurements  and  predictions  of  tangen¬ 
tial  velocity  in  the  swirl  chamber  with  the  probe  located  at 
10  mm  from  the  wall  surface  at  plane  2  which  is  45°  _o  the  port 
axis.  The  solid  line  represents  the  predicted  values.  A 
quite  sufficient  agreement  is  noticed  BTDC  as  shown  in  Fig. (30) 
however  discrepancy  is  noticed  during  the  mid  of  the  compres¬ 
sion  stroke,  Fig. (31).  This  discrepances  may  be  attributed  to 
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the  error  resulting  from  the  determination  of  the  inclination 
of  the  velocity  vector  by  the  hot  wire  anemometer.  The  measur¬ 
ed  values  are  clearly  sensitive  to  this  inclination  angle  which 
can  be  defined  accurately  after  the  circulating  flow  has  been 
established. 

It  is  planed  to  investigate  the  effect  of  various  geo¬ 
metry  parameters  of  the  swirl  chamber  and  working  conditions 
on  the  flow  field  to  isolate  some  of  the  important  features 
of  flow  and  to  validate  the  computational  procedure.  The 
present  work  lays  the  framework  for  a  more  detailed  study  in¬ 
cluding  spray  mixing  and  combustion  which  might  well  prove 
helpful  to  the  designer. 
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CONCLUSIONS 


The  multifuel  engine  must  operate  in  the  same  way  when 
using  either  clear  fuels  or  mixtures  with  an  arbitrary  ratio. 
An  analysis  of  the  influence  of  physical  properties  of  fuel 
on  mixture  creation,  combustion  and  engine  performance  has 
been  carried  out.  It  is  concluded  that  the  delivery  pressure 
must  be  increased  to  1.5  up  to  3  bar  and  an  electric  deli¬ 
very  pump  whicn  is  independent  on  the  drive  from  the  engine 
must  be  used  to  limit  the  vapour  lock.  The  injection  pump 
must  be  equipped  with  an  arrangement  for  changing  the  inject¬ 
ed  quantity  in  dependence  of  the  fuel  specific  weight.  Three 
types  of  combustion  chambers,  which  are  mainly  precombustion 
chambers,  are  recommended  to  use  in  multifuel  engines.  The 
analysis  leads  to  the  necessary  study  of  the  turbulence  mo¬ 
tion  of  air  in  combustion  chamber.-  heat  transfer  to  the  cham¬ 
ber  and  droplets,  droplet  size  distribution,  rate  of  fuel 
evaporation  and  the  rate  of  chemical  reaction.  The  study  of 
the  mentioned  parameters  are  required  to  satisfy  a  sufficient 
accuracy  of  the  spray  modelling. 

A  series  of  experiments  has  been  carried  out  on  a  spe¬ 
cial  designed  swirl  chamber  diesel  engine  to  investigate  the 
effect  of  swirl  chamber  geometry  and  operating  conditions  on 
the  instantaneous  surface  heat  transfer.  It  is  concluded  that 
the  mean  surface  temperature  decreases  with  the  increase  of 
engine  speed  and  decrease  of  load, swirl  chamber  diameter  and 
throat  diameter  of  the  tangential  port.  It  is  also  found  that 
the  instantaneous  heat  transfer  coefficient  to  the  wall  of  the 
swirl  chamber  is  described  quite  accurately  by  the  following 
dimensionless  equation; 

Nu  =  1.175xlO~4Re°*5Pra4  Eu_1 ' 5  (ds/dfc) 1 ' 2 

It  is  concluded  also  that  high  engine  speeds  and  greater  va¬ 
lues  of  swirl  chamber  diameter  results  in  high  heat  transfer 
to  walls. 

A  measuring  technique  and  a  prediction  procedure  for 
flow  in  motored  diesel  engines  swirl  chambers  are  presented. 
The  prediction  procedure  is  based  on  solving  numerically  the 
differential  conservation  equations  for  a  transient  turbulent 
flow  with  heat  transfer.  The  obtained  agreement  between  pre¬ 
dictions  and  measurements  is  fairly  good  which  confirms  the 
accuracy  of  the  prediction  procedure.  However,  the  prediction 
procedure  is  currently  being  tested  for  variables  more  rele¬ 
vant  to  turbulence  modelling  such  as  the  turbulence  intensify 
and  the  turbulent  shear  stresses. 
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